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FOREWORD 


This  report  was  prepared  by  the  Research  Division,  ManLabs,  Inc. , 
under  USAF  Contract  No,  AF33(657)-9826.  This  contract  was  initiated  under 
Project  No.  7350  "Refractory,  Inorganic  Non-Metallic  Materials",  Task 
735001  "Non-Graphitic" .  The  work  was  administered  under  the  direction  of  the 
A.F.  Materials  Laboratory,  Research  and  Technology  Division  with  Mr.  Fred 
Vahldiek  acting  as  project  engineer. 

This  report  covers  the  period  of  work  from  September  1962  to  September 

1963. 

Personnel  participating  in  the  work  included  L.  Kaufman,  H.  Bernstein, 
A.  Sarney  and  L.  Lindonen. 


ABSTRACT 


The  thermodynamics  of  interstitial  solid  solutions  and  non-stoichio  - 
metric  compounds  have  been  applied  to  analyses  of  the  Ti-C,  Zr-C,  Hf-C, 
Nb-C,  Ta-C,  Ti-O,  Zr-O,  Ti-N,  and  Zr-N  systems.  The  thermodynamic 
framework,  based  on  the  Schottky-Wagner  model  for  non-stoichiometric 
phases  has  been  applied  to  predict  phase  equilibria  and  vaporization  data 
for  these  systems.  Comparison  with  experimental  thermodynamic,  vaipor 
pressure,  and  equilibria  data  available  in  the  literature  is  very  satisfactory. 
In  addition,  correlation  between  thermodynamic  properties  and  vacancy  con¬ 
centration  data  where  possible  yields  good  agreement,  Results  of  x-ray 
measurements  of  the  expansion  coefficients  and  Debye  temperatures  of  ZrC 
and  HfC  performed  over  the  temperature  range  25°C  to  1500°C  are  also 
presented. 


This  technical  documentary  has  been  reviewed  and  is  approved. 


W.  G,  Ramke 

Chief,  Ceramics  and  Graphite  Branch 
Metals  and  Ceramics  Division 
AF  Materials  Laboratory 
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1.0  Introduction  and  Summary 

During  the  past  twenty  years,  increasing  emphasis  has  been  placed  on  the 
determination  of  the  thermochemical  properties  of  refractory  compounds.  As 
experimental  research  programs  designed  to  yield  information  on  the  integral 
free  energy  of  formation  of  these  compounds  developed,  the  experimental  dif¬ 
ficulties  characteristic  of  such  high  temperature  measurements  became  increas¬ 
ingly  apparent.  Moreover,  the  experience  gained  in  performing  reliable  meas¬ 
urements,  has  yielded  considerable  insight  into  the  costs  of  such  investigations. 
Under  these  circumstances,  it  is  quite  apparent  that  a  simple  theoretical  de¬ 
scription  of  the  thermodynamic  properties  of  refractory  compounds  which  could 
integrate  the  phase  equilibria,  vaporization,  and  integral  free  energies  of  com¬ 
pound  phases  would  be  extremely  useful.  The  development  of  such  a  theoretical 
model,  described  in  terms  of  explicit  numerically  predictable  parameters,  could 
be  useful  in  correlating  diverse  data  and  in  providing  sound  methods  for  estimating 
measurable  quantities  prior  to  the  more  expensive  experimental  studies. 

The  research  reported  in  the  body  of  this  report  describes  the  efforts  made 
along  these  lines  during  the  past  year.  The  theoretical  model  chosen  for  sim¬ 
plicity,  to  describe  the  properties  of  NaCl  type  refractory  carbide,  oxide,  and 
nitride  phases  was  the  Schottky-Wagner  model  which  was  first  proposed  about 
30  years  ago.  This  model  has  been  applied  for  the  specific  case  of  a  two-sub- 
lattice  compound  which  might  exist  over  a  wide  range  of  composition  and  tempera¬ 
ture.  By  exploring  the  means  of  applying  this  model  to  the  esse  at  hand,  it  has 
been  possible  to  integrate  phase  equilibria,  vaporization  information,  and  integral 
free  energies  of  formation  already  available  for  r efracto ry  compound  systems 

Manuscript  released  by  the  authors  8  October  1963  for  publication  as  an  ASD 
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in  addition  to  specifying  direct  numerical  methods  for  predicting  data  which 
is  presently  not  available  .  Although  the  present  model  is  simple,  possibly 
naive,  it  is  explicit  and  offers  a  method  for  direct  calculations. 

The  systems  which  have  been  considered  are  Ti-C,  Zr-C,  Hf-C,  Nb-C, 
Ta«C,  Ti-O,  Zr-O,  Ti-N  and  Zr-N.  Although  the  present  method  of  treating 
the  stability  of  refractory  compounds  is  not  represented  to  be  a  panacea  or 
a  substitute  for  further  experimental  studies,  its  usefulness  is  quite  apparent. 
Example  comparisons  of  experimental  and  computed  phase  diagrams,  vapor 
pressure  vs.  composition  curves,  congruent  vaporization  phenomena,  and 
thermodynamic  properties, detailed  in  this  report,  are  quite  encouraging. 

A  substantial  part  of  Section  2  of  this  report  will  be  published  as  a  portion 
of  a  paper  by  Larry  Kaufman  and  Edward  V.  Clougherty,  entitled  "Thermo¬ 
dynamic  Factors  Controlling  the  Stability  of  Solid  Phases  at  High  Temperatures 
and  Pressures".  This  paper,  which  was  presented  at  the  Winter  Meeting  of  tho 
American  Institute  of  Metallurgical  Engineers  in  Dallas,  Texas  on  26  February 
1963,  will  be  published  in  a  volume  entitled  "Metallurgy  at  High  Temperatures 
and  High  Pressure"  by  AIME  during  1964. 

In  addition  to  the  thermodynamic  investigation  described  above,  x-ray 
studies  of  thermal  expansion  and  lattice  vibrations  have  been  performed  on  ZrC 
and  HfC.  Measurements  of  Debye  temperatures  and  coefficients  of  expansion 
for  these  compounds  have  been  carried  out  at  temperatures  between  25°C  and 
1500°C  using  monochromatic  radiation.  The  purpose  of  this  study  is  to  permit 
intercomparison  between  x-ray  Debye  temperatures  and  Debye  temperatures 
derived  from  specific  heat.  These  compounds  were  chosen  for  study  because 
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of  their  practical  interest  and  for  the  additional  reason  that  while  reliable  low 
temperature  specific  heat  data  exist  for  ZrO  no  such  data  are  presently  avail¬ 
able  for  HfC. 

The  x-ray  results  given  in  section  3.0  were  used  to  estimate  the  vibration 
free  energy  of  HfC  and  its  contribution  to  the  free  energy  of  formation  of  HfC. 
The  estimated  free  energy  of  formation  of  HfC  is  given  as  a  function  of  tempera¬ 
ture  in  TABLE  IX.  This  latter  information  was  used  in  the  computation  of  phase 
equilibria  (p  22)  and  vaporization  (pp  29,  30,  36,  and  37)  in  the  hafnium -carbon 
system  performed  in  section  2,0, 
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2,  0  Thermodynamic  Factors  Controlling  the  Stability  of  Solid  Phases  at 
High  Temperatures 

Delineation  of  the  thermodynamic  factors  controlling  polymorphism  in 

pure  metals,  and  phase  equilibria  in  binary  substitutional  metallic  solutions 

(1-51 

has  been  accomplished  with  some  success'  The  purpose  of  this  factori¬ 

zation  is  to  predict  and  interpret  high  temperature  phase  equilibria  in  metallic 

:  (6-15) 

systems 

In  order  to  extend  this  approach  to  consideration  of  interstitial  solid  solu¬ 
tions  and  intermetallic  compounds,  it  is  necessary  to  perform  an  analogous  fac¬ 
torization  of  energy-entropy  terms  and  their  respective  concentration  dependences 
which  characterize  these  phase  classifications.  In  particular,  it  is  necessary  to 
recognize  the  interaction  between  phase  equilibria  and  free  energy^^,  and  to 
formulate  the  free  energy  of*icompound  phases  In  a  manner  which  implicitly 
presupposes  their  "non-stoichiometric"  nature^1^. 

Fortunately,  statistical  models  for  these  situations  have  been  presented  in 
detail  so  that  the  present  discussion  need  deal  only  with  the  application  of  these 
models  to  specific  cases.  Although  these  models  are  idealized,  they  offer  a  ra¬ 
tional  first  approximation  to  the  present  problem.  The  Schottky-Wagner^^'  * ^ 
model  of  non-stoichiometric  phases  will  be  used  as  the  framework  for  describing 

the  properties  of  compound  phases,  while  the  interstitial  solutions  will  be  de- 

(41 

scribed  in  terms  of  a  model  presented  earlier  '. 

As  a  final  statement  of  introduction  it  is  necessary  to  define  the  basis  of 
the  present  discussion  with  reference  to  the  size  of  the  system  and  the  reference 
states.  Following  the  suggestions  of  Wagner^ ^  and  Kubaschewski^^'  it 
seems  most  convenient  to  restrict  discussion  to  a  system  containing  Avogadro's 
number,  N,  of  atoms.  Hence  the  mass  basis  is  one  gram  atom  of  solution  or 
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one  gram  atom  of  compound  phase.  Accordingly  the  composition  of  a  given 
phase  is  specified  by  stating  the  atomic  fractions  of  the  component  elements. 
Although  this  system  is  by  no  means  universal  it  appears  (in  the  opinion  of 
the  present  authors)  to  be  most  convenient. 

(17) 

In  concurrence  with  the  viewpoint  of  Searcy'  ',  the  reference  states  are 

chosen  such  that  the  free  energies  of  the  pure  elements  in  their  stable  form  at 

one  atmosphere  pressure  and  0°K  are  zero.  Although  the  choice  of  such  a 

standard  state  is  subjective,  the  present  authors  heartily  endorse  the  suggestion 
(17) 

that'  '  "The  usual  thermodynamic  convention  of  assignment  of  an  activity  of 
unity  to  the  compound  cripples  efforts  to  understand  the  fundamental  behavior 

(19) 

of  the  system.  The  convention  of  referral  always  to  activities  of  the  thermo¬ 
dynamic  components  of  the  compound,  on  the  other  hand,  is  a  clear  and  uncom¬ 
plicated  framework  of  the  detailed  properties  of  the  compound". 

Since  the  symbolism  of  solution  thermodynamics  is  by  no  means  universal, 
a  Table  of  Symbols  is  appended. 

The  simplest  case  to  be  considered  along  the  lines  discussed  above  is  that 
of  a  binary  system  B-C  in  which  there  is  a  single  compound  existing  over  a  range 
of  composition  near  the  equiatomic  composition  (i.  e,  ,-Cq  g)  along  with  very 
little  solubility  of  C  in  B  or  B  in  C.  If  in  the  temperature  range  of  interest, 
the  compound  exhibits  the  (J  crystal  structure  while  B  and  C  represent  the  (3 
and  y  structures  respectively  then  the  B-C  phase  diagram  is  exemplified  by  the 
Ti-C  and  Zr-C  cases,  where  B  stands  for  Ti  and  Zr,  y  is  the  graphitic  form  of 
carbon,  C  is  the  sodium  chloride  lattice  arrangement,  and  x  is  the  atomic  frac¬ 
tion  of  carbon. 
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Figures  1  and  2  show  current  versions  of  the  Ti-C^^  and  Zr-C^^ 
phase  diagrams  as  well  as  "calculated"  or  idealized  phase  diagrams  for 
these  systems.  These  "calculated"  diagrams  (which  will  be  discussed  in 
detail  below)  are  not  presented  to  correct  or  to  cast  doubt  on  the  existing 
diagrams  but  rather  to  illustrate  in  a  semi-quantitative  way  the  relations 
between  the  free  energies  of  the  phases  involved  and  the  phase  stability. 

Figure  3  shows  the  integral  free  energy  of  ZrC  as  a  function  of  com¬ 
position  at  2000°K.  In  addition  Fig.  3  shows  the  partial  molar  free  energies 
of  zirconium  and  carbon,  which  are  obtained  graphically  by  tangent-intercept 
construction  or  mathematically  from  E^.  (Al-7).  Fig.  3  has  been  calculated  on 
the  basis  of  the  idealized  model  which  will  be  discussed  below.  The  lower  panel 
of  Fig.  3  shows  the  compositional  variation  of  the  partial  free  energies  of  zir¬ 
conium  and  carbon  across  the  ZrC  O’  field.  Note  the  large  variation  with  com¬ 
position.  If  the  pressure  of  zirconium  over  pure  /3Zr  (under  atmospheric  con¬ 
ditions)  is  p^^  and  the  pressure  of  zirconium  over  ZrC  of  composition  x  is* 

p2rtxl  thcn> 

rZr[xl  "  *Zr  =  RT  ln  PZr[x]/4r  (1) 

where  is  the  partial  molar  free  energy  of  zirconium  at  a  composition 

x  and  F^r  is  the  free  energy  of  (3Zr.  As  indicated  earlier  the  reference  state 

is  FZr[0°K]  =0.  Similarly! 

- - - - - 

In  this  paper,  square  brackets  will  be  used  to  denote  functional  relations.  Thus, 
PZr[x]  is  read  the  pressure  of  zirconium  over  the  O'  phase  which  is  a  function 
_of  composition. 

tFor  purposes  of  the  present  discussion  the  higher  polymers  of  gaseous  carbon 
i.  e.  Cg,C3,C^,  and  Cg  etc.,  will  not  be  considered  in  the  present  case  although 
they  are  certainly  important  in  the  temperature  range  to  be  considered.  Allow¬ 
ance  for  the  association  of  gaseous  carbon  atoms  can  be  made  by  utilization  of 
available  tabulations  of  thermodynamic  functions^^"^”)  for  the  data  of  the  free 
energies  of  association. 
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Fig.  I  — The  Zirconium-Carbon  phase  diagram. 
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Fig.  3  -  Calculated  integral  and  partial  free  energy  curves  for  the 
Zr-C  system  at  2000°K. 
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f£[x]-f£  =  RT  lnp£[x]/p£ 


(2) 


Reference  to  Fig.  3  shows  that  a  large  variation  of  p^J^x]  »  P^tx]  '  and 
the  ratio  P^/Pc  (dr°ppk*g  the  [*])  t0  he  expected  across  the  <J  field.  Of 
particular  interest  is  the  composition  of  the  free  energy  minimum  where 
9F°/9x  equals  zero.  At  this  composition  (Eq.  Al-7),  *  1^5.  It  should 

be  noted  that  this  composition  need  not  coincide  with  the  stoichiometric  com¬ 
position. 

Figure  4  shows  free  energy-composition  curves  for  the  (3,  y.  O'  and  L 
(liquid  phases)  of  the  Zr-C  system  at  temperatures  above  and  below  the  melting 
point  of  Zr-C  (which  is  assumed  to  occur  at  x  =  0.5)  and  above  and  below  the 
other  invariant  temperatures.  These  curves  are  idealized  since  they  are  based 
upon  several  simplifications  which  will  be  discussed  in  detail  later  but  they 
should  be  illustrative  of  the  general  features  of  the  real  case.  Fig.  4  shows  the 
construction  describing  the  free  energy  of  formation  of  stoichiometric  ZrC  i.  e. 
AF^[0.5]  or  AF^(dropping  the  [0.5])  in  cal/g.atom.  In  this  case  2AF^  is  the 
free  energy  of  formation  in  cal/mol  which  is  the  generally  tabulated  func- 
tion^^*^,  Fig.  4  illustrates  the  graphical  "common-tangent"  construction 
which  determines  the  phase  boundaries.  The  mathematical  equivalent  of  the 
common-tangent  construction  is  the  equilibration  of  the  partial  molar  free 
energies  of  both  components  across  the  two  phase  field.  For  example  at  T  = 
3400°K  in  the  ZrC  system  (Figs.  1  and  4)  the  partial  molar  free  energy  of 
zirconium  in  the  liquid  phase  at  a  composition  x^^-  is  equal  to  the  partial  molar 
free  energy  of  zirconium  in  the  <J  phase  at  a  composition  XgL  and  a  similar 
equation  holds  for  carbon  atoms.  Mathematically, 
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and 


(4) 


Eqs,  3  and  4  are  the  equivalent  of  the  rule  of  "common-tangents". 

One  of  the  general  features  of  Fig.  4  which  does  not  depend  upon  models 
or  approximations  are  the  relative  magnitudes  of  AF^  F^,  and  F^y  as  functions 
of  temperature.  Since  the  entropy  of  zirconium  is  much  larger  than  that  of  car¬ 
bon  at  any  temperature  it  has  a  much  lower  free  energy.  Coupling  this  fact  with 
the  known  compositional  limits  of  the  O’  phase  (i.  e,  x^s  0.5  and  Xg^  less  than 
0. 5)  require  that  the  composition  at  which  F  is  a  minimum  must  be  located  at 
x<  0,  5. 

In  order  to  inquire  more  deeply  into  the  factors  determining  the  tempera¬ 
ture  and  compositional  range  of  stability  of  compound  phases  it  is  necessary  to 
consider  the  "shape"  of  the  free  energy  vs.  composition  curve  (i.  e.  as  shown 
in  Fig,  3)  and  the  factors  which  determine  the  shape  in  the  vicinity  of  the  mini¬ 
mum.  Such  considerations  are  best  demonstrated  by  discussing  specific  exam¬ 
ples.  The  examples  which  will  be  discussed  in  subsequent  portions  of  this  section 
are  the 

a.  Ti-C,  Zr-C,  Hf-C,  Nb-C,  and  Ta-C  systems 

b.  the  Ti-O  and  Zr-O  systems 

c.  the  Ti-N  and  Zr-N  systems 

As  indicated  earlier,  the  consideration  of  these  systems  will  be  made  on 
on  the  basis  of  the  Schottky-Wagner  model  in  whose  context  these  cases  (a,  b,  and 
c)  represent  three  distinct  classes.  In  some  respects  however,  these  systems 

are  similar  since  all  the  compounds  of  interest  are  refractory  phases  of  the 

.  # 

sodium  chloride  type  which  exist  over  a  range  of  composition  near  x  =  1/2  . 

5“*” 

The  Zr-O  system  does  not  exhibit  the  sodium  chloride  x  =  1/2  structure. 
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Similar  analyses  can  be  applied  to  other  phase  classifications  having  different 
crystal  structures  corresponding  to  other  stoichiometric  compositions.  The 
difficulties  in  dealing  with  other  stoichiometries  (i.  e.  1/3, 1/4,  E/5  etc.)  are 
algebraic  rather  than  conceptual  (See  Appendix  Al). 

Although  a  detailed  description  of  the  particular  variant  of  the  Schottky- 
Wagner  model  used  in  the  present  considerations  is  given  in  Appendix  Al  it  is 
of  value  to  present  the  reader  with  a  qualitative,  quasi-physical  description  of 
the  model  before  launching  into  a  detailed  application  to  specific  systems.  The 
best  point  of  departure  for  such  a  consideration  is  the  vacancy  parameter  or 
which  is  directly  proportional  to  the  fraction  of  vacant  lattice  sites  at  stoi¬ 
chiometry.  For  the  case  of  50/50  compounds  2a  is  the  fraction  of  metal  (or 
non-metal)  lattice  points  which  are  vacant.  If  Qf  is  very  small  l.  e.  less  than 
0, 1%  then  (for  50/50  compounds) 

2AF°[  x  =0.  5]  =  -FMe+-  Fc+  -  2  RT  In  2  a  (5) 

where  AF^  is  the  free  energy  of  formation  of  the  stoichiometric  compound  in 
cal/g.  atom.  In  2q.(5),  FMq+  and  F^  are  the  free  energies  of  formation  of 
metal  and  carbon  vacancies  in  the  compound.  Consequently,  three  new  param¬ 
eters,  FMe+I  F^.+r  and  a  are  introduced  and  interrelated  through  Eq,  (5).  These 
parameters  control  the  "shape"  of  the  free  energy  curve  near  its  minimum  point. 
Thus,  on  the  metal  rich  side  of  stoichiometry  there  is  an  excess  of  C  vacancies 
and  the  rate  of  ascent  of  the  F^vs  x  curve  will  depend  on  the  magnitude  of  F^, 
On  the  carbon  rich  side  the  reverse  is  the  case.  Since  the  partial  free  energies 
in  a  sense  reflect  the  compositional  variation  of  F^[x]  they  are  determined  by 

bV  FMe+'  FC+’  and  a* 
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Finally,  Eq.{5)  states  that  for  a  given  free  energy  of  formation,  the  frac¬ 
tional  number  of  vacancies  will  depend  upon  the  free  energy  of  formation  of 
metal  and  carbon  vacancies.  This  equation  is  obtained  by  writing  a  general 
equation  for  the  free  energy  of  a  gram  atom  of  compound  of  a  fixed  composi¬ 
tion  at  a  fixed  temperature  and  pressure.  This  compound  is  envisioned  to  con¬ 
tain  an  arbitrary  number  of  metal  and  non-metal  vacancies  which  contribute 
positive  free  energy  increments  through  Fj^e+  anc*  FC+  and  negative  free  energy 
increments  via  the  entropy  of  mixing.  The  overall  free  energy  is  then  minimized 
(at  constant  temperature,  pressure  and  composition)  by  allowing  the  number  of 
vacancies  to  vary,  yielding  equation  (5)  as  the  criterion  for  minimum  free  energy. 

2.  1  Discussion  of  the  Ti-C,  Zr-C,  Hf-C,  Nb-C  and  Ta-C  Systems 

In  order  to  proceed  with  the  discussion  of  the  metal-carbon  systems,  it  . 
is  necessary  to  consider  the  properties  of  the  graphite  form  of  carbon,  C7 ,  which 
sublimes  at  about  4000°K^2^ .  By  combining  the  thermodynamic  data  for  C7  with 
the  entropy^20'*  and  enthalpy^2^  of  the  diamond  form  of  carbon,  C^,  at  (298 °K) 
and  the  high  temperature  data^2^>  2^  for  and  C7,  the  free  energy  difference 
between  6  and  y  can  be  approximated  (for  T>400°K)  as 

Af£"7  «  -280  -  1,1 6T  cal/g.at.  (6) 

Since  the  volume  change  at  room  temperature  can  be  computed  from 
the  volume  per  gram  atom  of  graphite^2®^  and  diamond^2^  as  AV^  7~  +  1.91 
cm2/g.  at.  The  pressure  dependence  of  Eq(6)  can  be  approximated  as^ 

AF^y[T,P]  2  -280  -  1.16T  +  23.9P(1.91)  cal/g.at.  (7) 

where  P  is  the  pressure  in  lcilobars*. 

Note  that  the  factor  23.9  in  Eq.  (7)  is  a  factor  converting  kilobars  to  cal/cm  . 
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Setting  Eq.  (7)  equal  to  zero  gives  a  reasonably  good  calculation 

(30  31\ 

of  the  T»P  curve  for  the  yZ 6  conversion  line  ’  . 

The  entropies  of  melting  of  the  diamond  cubic  forms  of  Si,  Ge,  and 
Sn^  are  6,6,  6.3,  and  5,3  cal/g.at.  °K  respectively.  Therefore  the  entropy 
of  melting  of  diamond  might  be  expected  to  be  about  6  cal/g.at.  °K.  Combina¬ 
tion  of  this  estimate  with  the  entropy  change  for  the  6— y  reaction  (Eq.  6)  results 
in  an  estimated  value  of  5  cal/g,  at  °K  for  the  entropy  of  melting  of  graphite. 
Bundy^3^  has  measured  the  enthalpy  of  melting  of  C^at  high  pressures  and 
temperatures  and  suggests  that  AS£,  =  5.4  cal/g.at.  Since  the  y-L-  vapor 
triple  point  occurs  at  about  0. 11  kbars  and  4000°K^3*'  3^  it  is  safe  to  estimate 
the  melting  point  of  graphite,  i.  e.  ,  as  about  4000°K.  Hence 

-  21,  600  - 5. 4T  cal/g.at.  (8) 

Although  this  expression  is  slightly  in  error  due  to  the  fact  that 
must  be  higher  than  the  sublimation  temperature  (since  graphite  sublimes 
at  atmospheric  pressure  and  4000°K)  examination  of  the  available  data^3*'  3^ 
indicates  that  the  difference  in  temperature  between  the  melting  point  and 
sublimation  point  is  not  large. 

In  order  to  calculate  the  Zr-C  and  Ti-C  phase  diagrams  the  following 
assumptions  are  made  (See  Figs  1  and  2} 

a.  x^g  and  x^L  *  0,  xy  g  and  xyL  =1,  and  xffy  and  xffL  *0.5 

b.  The  melting  point  of  the  Cl  phase  occurs  at  x  *  0,5 
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These  "Idealizations"  are  not  too  far  from  reality  and  are  used  to 
simplify  the  arithmetic.  Finally  a  more  questionable  assumption  is  made  con¬ 
cerning  the  properties  of  the  liquid  phase  i,  e,  that  it  is  regular.  Hence, 

FL  «<l-x)F^e  +  xFq  +  RT  (xlnx+(l-x)ln(l-x))+Lx(l-x)  (9) 

where  L  is  a  constant,  The  only  palatable  feature  of  this  assumption  is  that 
L  is  directly  specified  since  at  x*0,5  and  T  ■  T^,  FL*Fff.  Assuming  that  a 
io  a  small  number  (i.  e.  that  AFff[0,5,T]'  AF^  (see  Table  Al-I)) 

L  ■  4  (AFa[Tff]  +TaRln2  +  0.5AF^''’r[TCJ])  (10) 

Using  the  tabulated  values^2^  for  AF^[  T]  (in  cal/g.at.),  Eq,(7), 
and  the  published  melting  points  3750°K,  3520°K,  and  4070°K  one  obtains 
L« -36,080,  -44,620,  and -42,  780  cal/g.  at,  respectively  for  the  Zr-C,  Ti-C, 
and  Ta-C  cases. 

Next,  advantage  is  taken  of  the  fact  that  x^s  0.5,  and  xy(j  =  1. 
Since  equilibration  across  the  O'  +  y  field  requires  that 


Then  (Eq.  Al-10) 


*c|xyff  "clxuy 

Fc+  »  -RT  In  2a 


(ID 


(12) 


The  x^g  and  St^g  vs  T  curves  can  now  be  calculated.  On  the 
basis  of  Eqs.  9,  10,  and  (Al-7), 

FMe8^e  +  RTln(1-x)+L*2  <13> 
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and 

b  f£  +  RT  lnx+L(l~x)2 

Equilibrating  the  partials  across  the  L  +  O’  fields  for  T<T^  (and  using 

Eqs,  Al-11-14)  yields 

(14) 

2AF°  +  AFj'*y  «  RTln(l-xLa)xLQ.  +  L  ((l-xLff)2  +  x^g)  +  RT  In  (l-x^x^ 
and 

2AFa+ AF^,'*V  «  RT  In  (1  -*Lff)  *Lg  +  L  ((1  -XLQ  )2  +  x2ff)  (15) 

where  AF^  and  AF^  y  are  understood  to  be  temperature  dependent.  Since 
the  last  term  on  the  right  of  Eq.  14  is  small  enough  to  be  neglected  in  the  first 
approximation,  Eqs.  (14)  and  (15)  can  be  used  to  calculate  XgL  andXg^.  For 
temperatures  below  Tjj^e  ,  must  be  added  to  the  left  side  of  Eq.  (14). 

This  quantity  is  approximately  equal^2^  to  -2.3  -  T) , 

The  calculation  of  the  x^  L  vs  T  curve  is  made  by  equilibrating  the 
metal  partials  across  the  |3  +  L  field  hence, 

AF^-RTlnd^pJ+Lx2^  (16) 

Similarly  the  xT  vb  T  curve  results  from 
*-•  r 

AFC~r  s  RT  lnxLy  +  L^-xLy)2  (17) 

The  eutectics  T^.  and  T^,  result  from  the  intersection  of  the  x^ 
and  x^g  curves  and  the  x^  and  curves  since  at  Tg,  =  x^g  while 

at  Tg,  x^y  *  x^g  .  Assumption  a)  above  requires  the  occurrence  of  eutectics. 

If  this  assumption  is  relaxed  (complicating  the  arithmetic)  then  peritectics  can 
occur. 
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Up  to  this  point  no  assumptions  concerning  the  0  phase  have  been 
made  other  than  Xpy  =  0,5  which  leads  to  Eq.  (12),  In  order  to  complete  the  calcu¬ 
lation  of  the  diagram,  an  additional  "post  facto"  assumption  concerning  the  O’  phase 
is  required.  This  additional  assumption  is  made  by  setting 

RTln  a  •  2AHff[  0°K]  (18) 

where  AH^[0°K]  is  the  enthalpy  of  formation  (cal/g,  at. )  of  the  compound  at 
0°K.  This  assumption  fixes  Fc+  ,  (Eq.  12),  and  FMq+  ,  (Eq.  5), 

FMe+  "  ^  °0K]  -  RT  In  2-2AFff[T]  (19) 

For  the  Zr-C  case  at  2000°K, substitution  of  the  tabulated  values^2^ 
yields  FZji+  *»  79,750,  F^,+  *>41,350,  and  RTlna  ■  -44, 100  cal/g.  atom  re¬ 
spectively.  Thus,  a  is  a  very  small  number  (i.  e.  about  10“5)  and  FZr+  is 
much  larger  than  FC+  leading  to  the  asymmetry  in  the  F^  vs.  x  curves  shown 
in  Figs.  3  and  4.  With  this  assumption,  the  Xg^  and  xff  L  curves  can  be  cal¬ 
culated  by  equilibrating  the  partials  of  the  metal  atom  across  the  (3+0  and 
0  +  L  fields.  The  results  are 

RT  In  (1  ~2xg^)  /  (l-x^)  *>  2AH°[  0°K]  -2AF°[T]  +RTln2  (20) 

for  T<T£  yielding  Xg^  vs  T.  The  calculation  of  XgE  is  made  in  the  same 
manner  for  T>Tg  from  Eq.  (21) 

RTln(l-2xgL)/  (l-XgL)  »2AH°[  0°K]  -2AF°[T]  +RTln2 

+  RTln(l-xLg)  +  x2LgL  (21) 
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for  T>^e.  If  T^e>T  >  T£  ,  must  be  added  to  the  right  side  of  Eq. 

(21).  Since  x^g  vs.  T  is  known  from  Eq.  (14)  at  this  point,  then  Xgj  vs.  T 
is  defined.  These  Xg^[T]  values  can  be  resubstituted  into  Eq.(14)  to  reitera- 
tively  recalculate  xLg{  T] . 

Briefly  then,  Figs,  1-4  result  from  assumptions  a)  andb),  and 
the  ad  hoc  assumptions  of  regularity  for  the  liquid  phase  and  Eq.  (18)  de¬ 
scribing  the  vacancy  concentration  of  the  0  phase.  Needless  to  say  these  equa¬ 
tions  can  be  applied  to  calculate  any  of  the  other  metal-carbon  phase  diagrams 
which  are  suspected  to  have  similar  features  i.  e.  Hf-C,  Nb-C,  V-C,  etc.  This 
might  be  a  useful  first  approximation  in  cases  where  phase  diagram  data  are 
absent  or  sketchy.  In  addition  these  equations  can  be  used  to  calculate  the  acti¬ 
vity  and  vapor  pressure  (i.  e.  Fig.  3)  across  the  single  phase  d  field. 

Fig.  5  shows  the  calculated  Ta-C  phase  diagram,  with  the  cal¬ 
culation  performed  as  if  Ta^C  is  non-existent.  Consequently,  the  modification 
which  might  be  expected  if  Ta2C  is  included  is  shown  by  adding  the  full  lines. 
The  observed  Ta«C  diagram  shows  the  Ta2C  peritectlc  near  3700 °K.  At  this 
temperature  the  calculated  value  of  x^g  is  about  0.  33  so  that  FL[  0.  33]  should 
be  approximately  equal  to  the  free  energy  of  Ta2C  at  this  temperature.  Per¬ 
forming  this  calculation  yields  a  free  energy  of  formation  of  Ta2C  at  3700°K 
as  -28,250  cal/g.  at,  while  the  tabulated  value  is  -29,  100  +  2500, 

Apart  from  using  this  formalism  to  develop  semi-quantitative 
expressions  for  the  partial  free  energies,  the  representation  of  the  composition 
dependencies  of  the  vapor  pressure  can  be  used  in  some  kinetic  considerations  . 
For  example  the  rate  of  Langmuir  evaporation,  G,  is  given  by 

G  =  44.  4  p  mols/cm^sec  (22) 
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The  composition  at  which  congruent  vaporization  occurs  at  a  given  tem- 
perature  xc,  is  defined  by  the  condition 

<23> 

or 

Pctxc1/pMe[xc]  SXc(1-!tcrlMMe/2MC1/2  <24> 

Since  p^[  xc]  /p^0[  XC1  can  ke  calculated  as  a  function  of  composition, 
xc  vs.  T  can  also  be  calculated  (see  Fig.  3). 

Equations,  1,  2,  5,  12,  18,  19,  A1  -9  through  12,  of  Appendix  AI  and  the 
known  vapor  pressures  of  pure  carbon,  Ti,  Zr,  Hf,  Nb,  and  Ta  tabulated  by 
JANAF  can  be  employed  to  predict  the  vapor  pressures  of  metal  and  carbon 
over  the  refractory  carbides.  The  pertinent  results  are: 

RT  lnPMe^  °* 55  /Pmo  a  -FMe+“  RT  ln  2" 

RT  lnpj  {[  0.  5]  /p°  *  -Fc+-  RT  ln  2a 
for  x»0.5  l.  e.  at  the  stoichiometric  composition,  and 

RT  ln  p£je[  x]  /p^Q  =  -FMq+  +  RT  ln  ( 1  -2x)  4“ 1  ( 1  -xja  *2 
RTlnpJ[x]/p°  =  -Fc+  +  RT  lnx(l-2x)”* 

when  x  is  less  than  0.5.  These  equations  can  be  applied  to  the  refractory  NaCl 

type  monocarbide  phases  in  the  Ti-G,  Zr-C,  Hf-C,  Ta-C  and  Nb-C  systems, 

0.  , 

since  the  free  energy  of  formation  at  stoichiometry  AF  [0.5]  is  known.  Note 
that  these  are  the  same  equations  used  in  computing  the  binary  phase  diagrams 
shown  in  Figs,  1-5.  The  following  sources  were  used  for  2AF^[0.5];  TiC  and 
ZrC  (Elliot  and  Gleiser,  Thermochemistry  for  Steelmaking  Addison-Wesley 
Publishers,  Reading,  Mass.,  I960):  -32,300  +  1.32T  cal/rnol  for  NbC  and 
-35, 900  +  1. 2T  cal/mol  for  T aC  (both  from  W.  Worrel,  ScD  Thesis,  M.I.T.  , 

x - - 

(Submitted  for  publication  to  J.  Physical  Chemistry  by  W.  Worrell  and  J.  Chip- 
man,  Oct.  1963,  In  the  publication  version,  the  authors  modified  the  NbC  equa¬ 
tion  to  -31, 800  +  1 . 6T  and  the  TaC  equation  to  -35,600  +  1. 6T.  These  modifica¬ 
tions  will  produce  minor  changes  in  Figs.  12-15). 
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May'63):  and  2AF°[0.5]  for  HfC  as  given  in  TABLE  IX*. 

These  calculations  can  be  performed  explicitly  at  any  temperature  desired 
and  are  shown  in  Figs. ,  6  through  15  at  selected  temperatures,  These  curves 
depict,  in  explicit  fashion,  the  vapor  pressures  of  both  components  over  the 
single  phase  refractory  carbide  field. 

Since  the  rate  of  Langmuir  evaporation  (Eq,  22)  can  be  expressed  in  terms 
of  the  vapor  pressure,  these  equations  or  curves  can  be  used  to  predict  rates 
of  evaporation.  Crossing  of  the  component  pressure  vs.  composition  curves 
{which  occurs  for  ZrC  and  HfC)  indicates  congruency.  Application  of  Eqs.(23) 
and  (24)  yields 

2  RT  In  ( 1  -2xc)  =  RT  In  pj/p^ "  0  •  5  RT  ln  FC+  +  FMe+  +  2  RT  ln  «  ( 2 5) 

where  x  is  the  composition  for  congruent  vaporization  and  M„  and  Mv .  the 
c  O  MO 

atomic  weights  of  carbon  and  metal  respectively.  Equation  (25)  explicitly  defines 
xc.  Once  xc  is  calculated  (as  a  function  of  temperature),  then  P^et  xc*  T]  and 
p^[xc,  T]  which  are  the  pressures  of  metal  and  carbon  at  the  congruently  vapori¬ 
zing  composition  as  a  function  of  temperature  can  be  predicted.  These  curves, 
predicted  for  ZrC  and  HfC  are  compared  in  Figs.  16  to  19  with  measurements 
at  G.  E.  For  TiC.where  congruency  is  not  indicated  by  Figs,  6  and  7,  the  rate 
of  evaporation  of  carbon  saturated  TiC  {1.  e,  x=0.5)  has  been  calculated  and 
compared  with  measurements  in  Fig.  20.  Excellent  agreement  is  obtained.  For 
the  TaC  and  NbC  cases,  there  are  no  vapor  pressure  data  available  to  compare 
with  the  predicted  pressure  composition  curves.  However,  E,  K.  Storms  and 

N. £ Uco  r  i  an_ ti*_  s  •_ 6em •  U%0)  64  1471)  indicate  that  in  niobium  carbide  at 

*In  the  calculated  Hf-C  phase  diagram  (i.  e.  similar  to  Figs  1, 2  and  5  with 
Twf  =2495°K,  =  4170°K)  the  HfC-C  eutectic  is  calculated  at  3100°K  and 

78ar/oC  while  the  HfC-Hf  eutectic  is  calculated  at  2200°K  and  7a/oC.  At  this 
temperature  the  Hf  rich  boundary  of  HfC  is  computed  at  36a/oC. 


-  22  - 


3070°K,  congruent  vaporization  occurs  at  a  composition  of  41.5  atomic  percent 
carbon.  The  predicted  value  (i.e,  see  Fig.  13)  is  about  40%. 

In  addition  to  the  demonstrated  value  of  the  present  approach  in  dealing 
with  the  prediction  of  phase  equilibria  and  vaporization,  there  is  an  additional 
area  of  importance  to  which  it  can  be  applied.  Berkowitz  (ASD  Technical  Report 
#62-203  pt  II  January,  1963  "Kinetics  of  Oxidation  of  Refractory  Metals  and  Alloys 
at  1000°C-2000°C"  pp  26-68)  has  applied  the  Webb,  Norton,  and  Wagner  theory  of 
oxidation  to  a  consideration  of  the  oxidation  resistance  of  group  IVA,  VA,  and 
VIA  metal  carbides.  Dr.  Berkowitz  has  shown  that  this  theory  can  usefully  pre¬ 
dict  an  upper  limit  to  the  temperature  at  which  these  carbides  can  be  applied.  In 
order  to  apply  the  WNW  theory  (pp  27-29,  Eqs.  II— X  through  II-9  of  TR  #62-203), 
it  is  required  that  thermodynamic  data  on  the  activity  of  carbon  and  metal  in  the 

compound  at  the  composition  in  question  be  available.  The  actual  numerical  values 
0*  O' 

of  ac[x]  and  aMe[  x]  play  an  Important  role  in  controlling  the  vapor  pressure  of 
volatile  oxidation  products  whose  formation  limits  the  operating  temperature. 

Thus,  the  application  of  the  Schottky-Wagner  model  provides  an  idealized 
but  explicit  framework  for  integrating  the  equilibrium  and  vaporization  properties 
of  the  carbides.  The  method  should  not  be  construed  as  a  replacement  for  experi¬ 
mental  research,  but  rather  as  a  guide  for  performing  the  extremely  difficult  high 
temperature  equilibria  measurements.  Moreover,  if  suitably  extended  to  the  three 
component  case,  it  can  be  used  to  make  intelligent  estimates  of  ternary  equilibria 
and  vaporization. 

A  final  point  of  application  of  this  method  is  in  the  consideration  of  me¬ 
tastable  equilibrium.  To  illustrate  this  point,  consider  Fig.  5,  which  shows  the 
calculated  Ta-C  phase  diagram  ignoring  TagC.  The  predicted  eutectic  at  2500°K 
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-  Calculated  compositional  dependence  of  the  vapor  pressures 
of  Ti  and  C  over  the  TiCPphase  at  2200 °K. 


z.  Calculated  compositional  dependence  of  the  vapor  pressures 
of  Ti  and  C  over  the  TiC0^  phase  at  3000 °K. 
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Calculated  compositional  dependence  of  the  pressure 
of  Zr  and  C  over  the  ZxZa phase  at  3000  °K. 


Calculated  compositional  dependence  of  the  vapor  pressure 
of  Hf  and  C  over  the  Hf  phase  at  22Q0°K. 
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-Calculated  compositional  dependence  of  the  vapor 
pressure  of  Hf  and  C  over  the  Hf  <Za  phase 
at  3000° K. 


u  -  Calculated  compositional  dependence  of  the  vapor  pressure  of  Nb  and  C 
over  the  NbC0^  phase  at  2200 °K. 


(sdJdMdsoujjo)  9jnssdJd  601 
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Fig.  i3  -  Calculated  compositional  dependence  of  the  vapor 
pressures  of  Nb  and  C  over  the  NbC0^  phase  at 
3000  °K. 


Calculated  compositional  dependence  of  the  vapor  pressures 
of  Ta  and  C  over  the  TaCa  phase  at  2200  °K. 


(S9JdqdS0tUJD)  9^nss9J<J  6o| 
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Calculated  compositional  dependence  of  the  vapor  pressures 
of  Ta  and  C  over  the  TaCa  phase  at  3000  °K. 


6  -  Comparison  of  calculated  and  observed  vapor  pressure 
of  Zirconium  over  congruently  vaporizing  ZrC. 


Comparison  of  calculated  and  observed  vapor  pressure 
of  Carbon  over  congruently  vaporizing  ZrC. 
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Comparison  of  calculated  and  observed  vapor  pressure  of 
Hafnium  over  congruently  vaporizing  Hf-C. 
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Comparison  of  calculated  and  observed  vapor  pressure 
of  Carbon  over  congruently  vaporizing  HfC. 
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(l/T°K)  xIO 

Fig,  2o  -  Comparison  of  calculated  and  observed  rates  of  evaporation  of  Titanium 
Carbide. 


ATOMIC- PERCENT 
OXYGEN 


Fig, 21  -  Lattice  parameters  and  phase  relations  in  the  titanium -oxygen 
(22,29,33) 
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system 


%o 


Fig.22 —  The  Zirconium- Oxygen  phase  diagram 
(after  Hansen  and  Anderko22). 
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Zr02  (Monoclinic) 


between  Ta  and  TaC  is  calculated  on  the  basis  that  TajC  does  not  form.  Recently, 
R.  V.  Sava  of  the  Union  Carbide  Research  Center,  Parma,  Ohio,  performed  an 
experiment  in  which  a  Ta-TaC  couple  was  heated  rapidly.  Ta,C  did  not  form  dur¬ 
ing  the  experiment  and  melting  was  observed  at  2200°C  (2470°K)  corroborating  the 
prediction  of  Fig.  5. 

2. 2  Discussion  of  the  Ti-O  and  Zr-O  Systems 

Figs.  21  and  22  show  the  Ti-0  and  Zr-O  phase  diagrams.  These 
differ  appreciably  from  the  Ti-C  and  Zr-C  situations  in  that  considerable  solu¬ 
bility  of  oxygen  in  both  the  Qf(h.  c.  p. )  and  (3(b.  c.  c. )  forms  of  these  metals  is 
exhibited  at  one  atmosphere,  several  "compounds"  are  formed  (in  Ti-O)  and  the 
CT  structure  (which  appears  only  in  the  Tl-O  system)  does  not  exhibit  the  high 
temperature  stability  indicated  in  the  Ti-C  and  Zr-C  cases;  moreover  the  most 
refractory  compounds  occur  at  values  of  x(atomic  fraction  oxygen)  which  are  dif¬ 
ferent  from  1  /2. 

Since  the  thermodynamic  properties  of  the  Ti-O  system  have  been 
(34-421 

most  extensively  evaluated'  ,  it  is  instructive  to  start  with  a  discussion  of 

this  system.  In  line  with  the  stated  reference  state  convention,  the  free  energy 

of  oxygen  (gaseous)  at  temperatures  above  500°K  and  one  atmosphere  will  be 
>!< 

approximated  by  Eq,  (25), 

j  400  -  4.  22  T  in  T/2.  74  cal/g.at.  (25) 

Kubaschewski  and  Dench^'*’  ^  have  measured  the  partial  molar 
free  energies  of  oxygen  for  several  Ti-O  alloys  in  the  interstitial  ^  and  O'  range. 
They  combined  these  data  with  published  values  for  the  free  energy  of  formation 
of  TiO^,  T^Oj,  and  higher  oxides  to  draw  curves  for  the  partial  molar  free 

T;  *“•“  “  . 

This  expression  differs  by  about  850  cal/g.at,  from  the  usually  tabulated  data'  ' 
which  gives  the  room  temperature  enthalpy  of  a  mol  of  oxygen  as  2075  cal/mol. 
This  value  refers  to  an  ideal  gas  between  0°K  and  ZOS^K.  In  spite  of  this,  the 
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energies  of  titanium  (F!^[x]  -F^j)  and  oxygen  (T^j  [  x]  -  1/2  Fq  )  at  1273°K, 
These  curves  are  naturally  consistent  with  the  tabulated  values  for  the  other 
compounds.  Mah,  Kelley  and  co-workers^^)  measured  the  enthalpy  of  forma¬ 
tion  of  a  Ti-O  alloys  at  298°K  for  x<  0.  25  while  Ariya,  Morozova,  and  Volf^37^ 
did  similar  measurements  over  the  range  0.20<x<0.70,  Mah  etal^^  also 

measured  the  specific  heat  of  two  alloys  containing  0.058  and  0.25  atomic  frac- 

(381 

tions  of  oxygen,  Hepworth  and  Schuhmann'  1  measured  the  pressure  of  hydro¬ 
gen  over  Ti-O-H  alloys  from  which  they  deduced  the  activity  of  oxygen  in  Ti-O 
alloys  at  1073°K  over  the  range  0<x<  0. 14,  These  calculations  were  made  on 
the  basis  of  an  infinitely  dilute  standard  state  and  then  converted  to  absolute 
activities  by  using  the  data  of  Mah  et  al^3^, 

(391 

Veinbachs,  Silver  and  Komarek'  '  measured  the  activity  of 

oxygen  in  the  Ti-O  system  over  the  temperature  range  1 100°K<T<1 300°K  for 

alloys  containing  up  to  25  a/o  oxygen.  They  used  their  data  to  compute  partial 

entropies  and  enthalpies  for  oxygen.  The  enthalpy  data  yielded  heat  of  solution 

in  the  dilute  solution  range  which  agree  with  the  calorimetric  values  of  Mah  et 

al.  In  this  range,  the  activity  data  agreed  with  those  of  Kubaschewski^3^’  33^and 
(381 

Hepworth'  .  However,  at  higher  compositions  the  results  of  Veinbachs  et 
(39) 

al'  differ  from  those  deduced  by  Kubaschewski,  Consequently  these  data  do 
not  agree  well  with  the  published  data  for  titanium  compounds. 


tables  give  the  298  K  entropy  (49.00  cal/mol)  based  on  solid  oxygen  at  0°K.  The 
enthalpy  at  298  K  based  on  solid  oxygen  at  0°K  is  nearer!^)  3750  cal/mol  which 
is  the  value  used  in  the  present  convention.  The  difference  between  3750  and 
2075  (divided  by  2  i.  e,  per  gm,  atom)  accounts  for  the  abovementioned  difference. 
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(39) 

On  the  basis  of  their  entropy  data,  Veinbachs  et  ai  were 
able  to  show  that  the  positional  entropy  of  a  Ti-O  alloys  was  consistent  with 
an  ideal  interstitial  solid  solution  (i.  e.  Eq.  A2-7,  n  =1  since  O'  is  h.c.p. 
and  z  =  1). 

Eqs,  (A2-1-8)  indicate  that  the  thermodynamic  properties  of  an 
interstitial  solid  solution  can  be  specified  if  several  parameters  are  known*. 
Thus,  HTio>  A,  and  0,pj[x]  need  to  be  specified.  The  latter  function  can  be  esti¬ 
mated  from  a  knowledge  of  the  atomic  volumes,  the  nr  citing  points,  and  the 
atomic  weights  of  Ti  and  O  (Fig.  6).  In  this  estimation  procedure,  is  taken^ 
as  1730°K,  The  resultant  0®^[x]  curve  is  approximately 

fijjtx]-  365  (l-x)(l+2x)  (26) 

This  equation  can  be  used  (with  A2-2)  to  estimate  the  specific 
heat  and  entropy  on  the  basis  of  the  two- Debye  0  method^  ignoring  contributions 
of  the  electronic  specific  heat&.  Comparison  of  the  calculations  with  the  reported 
specific  heat  data  50°K<  T<  300°K(36>  is  within  3%  for  x  =  0.058  and  10%  for 
x  =  0.25.  The  room  temperature  entropy  values  for  x  =  0,058  are  6,78  cal/g,  at. 
°K  (observed)  and  6.6  (calculated)  while  for  x  ■  0.25,  5.  15  is  observed  and  5.  50 
is  calculated.  Although  Mah  et  al^3^  suggest  that  the  excess  vibrational  entropy 
is  zero  (i.  e.  the  partial  vibrational  entropies  of  Ti  and  O  are  independent  of 
composition)  the  abovementioned  estimation  procedure  yields  concentration  de¬ 
pendent  partials.  This  arises  from  the  pronounced  maximum  in  the  melting 
point  at  x  =  0.  25,  It  would  be  of  interest  to  have  measurements  of  the  specific 


^Both  z  and  n  are  unity  as  indicated  above. 

%  ,  (40)  o 

Low  temperature  specific  heat  measurements'  1  on  a  Ti-O  alloys  (1-4  K) 

show  the  electronic  specific  heat  to  be  small. 
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heat  (50-300°K)  at  a  higher  oxygen  content  within  the  a  field  (0.25<  x  <  0.33) 
in  order  to  resolve  this  point. 

The  remaining  parameters  H^.q[0°K]  and  A  can  be  directly 
evaluated  from  Kubashewski*s  values  for  the  compositional  dependence  of 
(Fq  -1/2F°  )  at  1273°K.  The  results,  A  =  -24,650  and  H^iO[0°K]  a  -60,000 
cal/g.at.  reproduce  the  published  results  to  within  1%.  The  latter  parameter 
leads  to  an  enthalpy  of  formation  which  closely  approximates  the  value  sug¬ 
gested^^  for  tjje  stable  form  of  TiO.  This  is  not  surprising  in  view  of  the 
structural  similarity  between  the  0  structure  and  a  h.  c.p.  metal  lattice  in 
which  all  of  the  interstitial  sites  are  filled.  Moreover,  extrapolation  of  the  a 
lattice  parameters  (Fig.  6)  to  x  ■  0.5  leads  to  a  volume  of  5.6  cm  /g.at.  vs 

3 

5.5  cm  /g.at,  for  the  0  structure.  On  this  basis  (ref.  20,  Fig.  93)  one  would 

expect  very  similar  heats  of  formation.  Finally  substitution  of  these  parameters 

into  Eq.  (A2-8)  permits  calculation  of  the  enthalpy  of  formation  of  a  Ti-O  alloys 

(3  3) 

which  agree  to  within  3%  (or  better)  with  the  values  measured  by  Mah  et  al  '. 
The  final  representation  of  the  Ti  and  O  partials  are  given  by  Eqs.  (27)  and  (28), 

F^.-F®.  »  RTln(l-2x)(l-x)-1  -49,  300x2  +  3  RT  (ln(l+2x)(l  -x)  (27) 

-x(l-4x)(l-x)‘1(l  +  2x)"1) 

and 

^-1/2Fq  * -1/2Fq  -120,000  +  3  RT  In  632/T  +  RT  lnx(l -2x)_1  (28) 

-24,  650(1 -4x+ 2x2)  +  3  RT(ln(l +2x)(l-x)+(l-4x)(l +2x)-1) 

As  a  final  check  on  Eq,  (27),  this  equation  can  be  used  to  calcu¬ 
late  vs,  T  for  comparison  with  Fig,  6,  This  calculation  is  performed  by 
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equilibrating  the  partials  of  Ti  across  the  «  +  0  field  and  applying  Eq,  A2-6 
to  the  /3  phase  (na3)  for  dilute  solutions.  The  AF,p.  H[  T]  values  used  for  this 
calculation  were  taken  from  reference  (1).  The  agreement  between  the  cal¬ 
culated  xap  vs.  T  curve  agrees  to  within  2  a/o  with  the  observed  curve  Fig. 

(6)  for  1300°K<  T  <  2000°K. 

Turning  to  the  TiO  O'  phase,  we  make  the  first  approximation  (as 
in  the  case  of  TiC,  ZrC,  and  TaC)  that  0/  is  very  small.  Consequently  Eqs. 
(Al-11-14)  can  be  applied.  Designating  the  TigOj  phase  by  the  symbol  the 
following  relations  can  be  written  at  1273°K. 


i  j  y.Q  z  *TTllxff(y  FTl"FTi+ +  RTln^1"2xff0f)//'la  1 

*o  >/2Fa,-  Fo+  +  RT  “  %<,/(■ -2xcra) 

where  =  0.  337  and  XgQ  =  0.  465  from  Fig.  6,  and 


(29) 

(30) 


tRTta(l-*ffxl/(2xn-l)  (31) 

1/2Fo2-Fot  +  RT1"l2x<n-1l/4“2l‘i)x  <«) 

where  XCT\~  0  .  56  and  x^g  is  near  0.  6. 

In  these  equations  and  are  given  by  Eqs.  27  and  28  while 
Kj.j+ »  Fq+  and  a  are  the  free  energies  of  formation  of  titanium  and  oxygen 
vacancies  and  the  vacancy  parameter  for  the  Ti-O  O’  phase  at  1273°K.  Since 
[  0.  337]  -  j  Fq  *  -90,  500  cal/g.  at.  at  1 273°K  from  Eq.  (28),  Eq.  30  yields 

F0+=  +95,  300  cal/g.at.  (33) 
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Multiplication  of  Eq,  31  by  0.  4  and  Eq.  32  by  0,  6  and  addition  yields 
AFX  [0,6,  1273]  a  0.8AFff[0.5,  1273]  -  0.  2F0+ -0.  4  RT  In  2or 

+  0.2RTln(2x(JK-l)(l-xffK)2x^  (34) 

Where  AFX  and  AF*^  are  the  free  energies  of  formation  of  the  "Ti^Oj"  and 
"TiO"  stoichiometric  phases  (cal/g.at.),  Substitution  of  the  proper  numerical 
values^2^  yields 

RT  In  2a  -  -  6500  cal/g,  at.  or  a -  0.05  (35) 

Thus  a  is  not  small  and  the  approximations  introduced  by  Eqs, 
29-32  are  not  valid.  These  equations  are  now  rewritten  using  the  more  accurate 
expressions  in  Table  Al-1  and  applied  over  the  range  1273°K<  T<  2000°K  to 
solve  for  a}F^,^+  and  Fq+  (since  we  can  write  four  equations  at  each  tempera¬ 
ture).  This  procedure  yields  the  following  results  for  the  TiO^phase 

FT.+  - 18, 100 -3.3  T  cal/g.at.  (36) 

Fq+-  119, 800  -19. 9T  cal/g.at.  (37) 

RT  In  a  /(I -2a)  ~  -8200  cal/g.at.  (38) 

On  the  basis  of  these  results,  we  can  approximate  the  enthalpy  of 

formation  of  Ti  and  O  vacancies  as  -  18.1  kcal/g.  at.  and  -  119.8  kcal/ 

g.at.  Substitution  of  these  results  into  Eqs.  (A1 -25-26)  permits  calculation  of 

the  compositional  dependence  of  the  enthalpy  of  formation  of  TiO^  for  compari- 

(37) 

son  with  the  results  of  Ariya  et  al'  ' ,  Comparison  of  the  calculated  and  ob¬ 
served  values  shown  in  Table  I  is  satisfactory. 

As  noted  previously,  the  value  for  the  vacancy  parameter  a  for 
the  TiO^  is  surprisingly  large  (Eqs.,  35,  38).  However,  this  is  quite  consistent 
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TABLE  I 


COMPARISON  OF  THE  CALCULATED  AND  OBSERVED*37* 
COMPOSITIONAL  DEPENDENCE  OF  THE  ENTHALPY  OF 
FORMATION  OF  TiO0^ 


X 

at,  fr. 

oxygen 

-AH  [ 
kcal/j 

obs.*37* 

298] 

at. 

calc. 

X 

at,  fr. 

oxygen 

-AHU  [  298] 
kcal/g,  at. 

obs*37* 

calc. 

0.447 

57.29*0.  14 

56.8 

0,502 

62.88*0.02 

63.0 

0.466 

60. 18*  0. 12 

59.0 

0.506 

63.42*0.02 

63.3 

0.484 

61.40*0.05 

61.0 

0.528 

65.43*0.02 

65.3 

0.492 

62.26*0.03 

61.9 

0.548 

67.34*0.04 

67. 1 

0.500 

(Interpolated) 

62.80 

62.8 

0.556 

0.560 

68.43*0.03 

68.86*0.02 

67.8 

68.4 

TABLE  II 

COMPARISON  OF  CALCULATED  AND  OBSERVED*41,  42)  VACANCY 
CONCENTRATIONS  IN  THE  TiO17  PHASE 


X 

Percentage  of  Vacant 

Percentage  of  Vacant 

at.  fr. 

Ti  Sites 

O 

Sites 

oxygen 

obs,*41'  42*calc.  (0(  *0.079) 

obs.*41,42> 

calc,  (a  *  0.079) 

0.41 

4 

7.0 

34 

35.2 

0.50 

15 

15.8 

15 

15.8 

0.53 

19 

21.4 

9 

11.4 

0.56 

23 

28.2 

4 

8.8 

0,57 

26 

30.6 

2 

8.0 

CALCULATED  TEMPERATURE  DEPENDENCE  OF  a  (Eq.  38) 

T°K  1273  1400  1600  1800  2000 

a  0.036  0.047  0.066  0.084  0.101 
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Increasing  Free  Energy  (kilocal  /gm.  atom) 


Atomic  Percent  Oxygen 


0 


f 


Fig,  23 


Note:  The  shapes, 
of  the  T^Oj.TijOJ 

Ti  02  free  energy 

curves  are  esti¬ 
mated, _ 
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Ti 
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Calculated  integral  and  partial  free  energies  in  the  Ti-0  system 
at  1600 °K 
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with  the  results  of  Ehrlich^41' 4^,  Ehrlich  determined  the  fraction  of  vacant 
Ti  and  O  sites  by  performing  lattice  parameter  and  density  measurements  on 
a  series  of  TiO^  alloys  quenched  from  high  temperatures.  His  results,  along 
with  calculated  vacancy  concentrations  (Eqs.,  38,  Al-1-3,  and  Table  Al-I) 
corresponding  to  1500°C  (a  a 0,079)  are  given  in  Table  II.  The  agreement  is 
quite  good. 

This  comparison  is  particularly  satisfying  since  it  offers  support 
for  the  physical  basis  of  the  model  since  the  calculation  of  a  is  performed  on 
thermodynamic  grounds  and  yields  values  of  a  which  compare  quite  favorably 
with  those  derived  from  physical  measurements. 

As  a  final  illustration  of  the  properties  of  the  TiO^  phase  it  is 
instructive  to  plot  the  integral  and  partial  free  energies  at  1600°K  as  a  function 
of  composition  using  Eqs,  (27),  (28),  (36-38),  and  the  equations  in  Table  Al-l, 
Values  for  the  free  energy  of  formation  used  throughout  are  those  tabulated  by 
Gleiser  and  Elliot 

It  is  interesting  to  compare  Figs.  3  and  4  for  ZrC^  with  Fig.  23 
for  TiO0^  in  order  to  examine  the  effects  of  the  parameter  Ot  and  the  free  energies 

of  formation  of  metal  and  nonmetal  vacancies  on  the  shape  of  the  F°-x  curve. 

»5 

Examination  of  these  figures  shows  that  a  smaller  value  of  a  (i.  e,  10  for  ZrC 

vs  10"1  for  TiO)  leads  to  a  "sharper"  F°-x  curve  and  a  larger  variation  in 

partial  free  energies  across  the  single  phase  field.  Moreover  when  the  free 

energy  of  formation  of  vacancies  of  one  of  the  species  is  small  relative  to  the 

O' 

free  energy  of  formation  at  stoichiometry  (i.  e,  Ti  in  TiO)  the  F  -x  curve  loses 
its  characteristic  minimum. 

The  Zr-O  system  (Fig.22)is  very  similar  to  Ti-O  system.  The 
chief  differences  are  the  slightly  lower  values  for  the  compositional  limit  of  the 
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a  phase,  the  higher  melting  point  of  ZrC>2  {relative  to  Tif^)  and  the  absence 

of  stable  compounds  between  the  a  phase  and  ZrC>2  and  ZrO^  (we  will  use 

these  phase  designations  for  the  tetragonal  and  monoclinic  forms  of  ZrC^). 

(44  29)  (T 

There  have  been  references1  ’  '  to  a  metastable  ZrO  .  All  of  these  ob¬ 

servations  are  interrelated  and  stem  from  the  greater  stability  of  Zr02  relative 
to  TiOj.  If  the  tabulated  free  energy  of  formation  of  ZrC^  were  used  in  place  of 
that  for  TiOj  in  construction  of  Fig,  23,  then  the  x  =  0,667  minimum  would  permit  a 
common  tangent  to  be  drawn  between  the  T  and  a  phase  below  the  CT,  T^O^ 

and  TijOg  phase  thus  leading  to  phase  relations  similar  to  those  shown  in  Fig. 22. 

(391  (451 

Veinbachs  et  al  and  Kubachewski  and  Dench'  1  have  measured 

the  activity  of  oxygen  in  the  a  Zr-O  phase.  Veinbachs  et  al  suggested  that  the 

positional  entropy  of  a  Zr-O  alloys  is  ideal.  Using  this  information  and  the 

procedure  described  earlier  to  estimate  the  properties  of  the  interstitial  Ti-O 

a  solid  solutions  yield  the  following  results: 

e2rU]  *  ~  260°K(1>,  0“[x]  *  621°K 

A  <*  -23,  800  cal/ g.  at,  and  0°K]  "  -63,  000  cal/g.  at. 

Comparison  of  these  values  with  the  Ti-O  parameters  discussed  earlier  shows  the 
similarity  between  these  systems.  The  analogues  of  Eqs.  (27)  and  (28)  for  the 
Zr-O  system  are, 

*Z^"FZr  =  &Tln(l-2x){l-x)_1  -47,600x2  cal/g.  at.  (39) 

ig-  1/2Fq  a  -1/2  Fq  -126,000  +  3  RT  In  621 /T 
2  2 

+  RT  In x{  1  -2x)  ” 1  -23,  800(1 -4x  +  2x2)  cal/g.  at.  (40) 
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These  equations  can  be  used  to  calculate  by  using  the 

tabulated  AFZr  ^  [  T]  values^.  As  was  the  case  for  the  Ti-0  system,  the 
vs  T  curve  is  calculated  in  agreement  with  the  published  vs  T  curve 
(Fig,  7).  Equilibrating  the  partials  of  Zr  and  O  across  the  Of  +  T  fields  yields  equa¬ 
tions  analogous  to  Eqs,  (29)  and  (30).  Assuming  Of  (the  vacancy  parameter),  for 
Zr02  is  small  yields  (1300°K<  T<2173°K) 

*Zrl^T  ■IZrl«Ta*  IZ,-SZrt  +*Tl»(2.3xTO)2/27a3(l-^r[I)2  (41) 

*ol^,T  "  '*  %-*(»*  RT1»«To/s-3*ra  <«) 

where  3AFT:  -  FZr+  -  2FQf-3RTln  3CX/41/3  (43) 

'J 

In  these  equations,  AF  is  the  free  energy  of  formation  of  stoi- 

j 

chiometric  ZrO ^  ,  F-  +  and  are  the  free  energies  of  formation  of  Zr  and 
T  1  /3 

O  vacancies  in  ZrC>2  ,  and  30f/4  '  is  the  fraction  of  vacant  sites  in  the  T  phase 
at  x  a  2/3. 

Eq.  (41)  can  be  solved  for  ^+[  T]  by  using  Eq.  40  and  the  phase 
diagram  yielding: 

FQ+=  128,400  -  17.  5T  (44) 

Eqs.  (39)  and  (41)  can  be  used  together  with  the  phase  diagram  to 
solve  for  the  sum  of  FZr+  +  3RTlna  and  the  result  checks  with  Eq.  (43). 

FZr+  +  3RTlna  =+3000  "l2-2T  (45) 

If  the  solubility  of  oxygen  in  ZrC>2(i.  e,  the  high  oxygen  phase 
boundary  of  the  Zr02  phase)  were  known,  then  a  (and  F^^)  could  be  estimated 
explicitly.  Ackerman,  Thorn  and  Winslow^4^  have  measured  the  vapor  pressure 
of  various  gaseous  species  over  ZrC>2  at  2750°K.  They  suggest  a  value  of  0?  near 
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pig.  24  -  The  Zirconium— Nitrogen  and  Titanium- Nitrogen  phase  diagrams 
(after  Hansen  and  Anderko22). 


-3 

10  at  this  temperature.  Hence  3RTlntt*  -113,100  cal/g.at.  and  FZji+5  116,100 
-  12. 2T.  This  result  implies  that  for  ZvO^  the  free  energies  of  vacancy  formation 
of  both  species  are  of  similar  magnitudes. 

2.  3  Discussion  of  the  Ti-N  and  Zr-N  Systems 

The  titanium -nitrogen  and  Zr-N  systems  are  quite  similar  to  each  other 
and  to  the  Ti-O  and  Zr-0  systems  in  some  respects.  Fig,  24  shows  the  established 
phase  diagrams  which  indicate  stabilisation  of  the  a  interstitial  solid  solutions  and 
refractory  O'  phases.  Since  there  do  not  appear  to  have  been  any  complete  studies 
of  the  thermodynamic  properties  of  the  Of  solutions,  the  thermodynamic  properties 
will  be  estimated  in  a  manner  similar  to  the  procedure  used  for  the  Ti-O  and  Zr-O 
a  solutions.  In  the  case  of  Ti-N  we  find  that  the  compositional  dependence  of  the 
Debye  temperatures  is  given  by  Eq.  (46). 


0“t[x]  «  365  (1+x),  e“[x]  =  675  (1+x)  (46) 

The  values  of  A  and  H^1N[0°K]  are  found  to  be  -18,600  and  -39,000 
cal/g.at.  respectively,  thus  the  partial  free  energies  of  Ti  and  N  are  approxi¬ 
mately  (l.  e,  Eqs.  A2-6-7) 

TTi"FTl  "  RTlnd-ExHl-x)-1  -37,  200x2  +  3  RT  (ln(l+x)-x(l+x)’1)  (47) 

and* 

^*-1/2f£  = -1/2F^  -78,000+  3  RTln675/T  +  RTlnx(l-2x)*1 
2  2 

-18,600(l-4x+2x2)  +  3  RT(ln(l+x)  +  (1 -x)(l+x) _1 )  (48) 


the 


xap 


As  in  the  case  of  Ti-O  and  Zr-O,  Eq.  (47)  can  be  used  to  calculate 
vs  T  curve  for  comparison  with  Fig.  24.  The  differences  are  less  than 


As  in  the  case  of  oxygen  (see  section  2.2)  the  values  used  above  are  slightly 
different  from  the  tabulated  values'2°)  due  to  the  choice  of  standard  states. 
This  difference  is  +(3300-2072)/2=614  cal/g.at. 
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0" 

2  a/o  nitrogen.  Assuming  the  a  for  the  TiN  phase  is  small  and  equilibrating 
the  partials  of  Ti  and  N  across  the  a  4  0  field  (i.e.  see  Eqs.  29,  30)  yields 


F..,  «  82,200  -  23.  3T 
N+  ’ 


(49) 


E^,i+  4  2  RTln  la  »  -3800  +  1.75T 

At  the  high  nitrogen  boundary  of  the  0  field,  Xq^,  =  l/^F^* 
hence  (A1  -14) 

Fn+  n -2  RTln  2a  4  RT  In  (2xffN-l)/(l-xff  N) 

Eqs.  (49)  and  (51)  yield 


(50) 


and 


RT  In  a  b  -24,700 


FT.+  s  45,  600  -  1.0T 


The  Zr-N  system  can  be  treated  in  a  similar  fashion  yielding 
02r[x]-265,  0^[x]~  665,  A  -  0  and  H^JO]  --42,000  cal/g.  at.  thus, 


(51) 

(52) 

(53) 


N 1 


ZrNL 


*Zr~FZr  ‘  RTln(1-2*H1-x)‘1  (54) 

Y®  -  1/2F^  *.  -l/2F°  -84,0004RTln665/T4RTlnx(l-2x)_1  (55) 

Calculation  of  the  xQ^  vs  T  curve  from  Eq.  54  agrees  very  well 

(47) 

with  Fig.  24,  Kibler,  Lyon,  and  DeSantis  have  measured  the  vapor  pressure 
of  nitrogen  over  the  ZrN^phase  at  1800<  T  <  2300°K  and  0.  415  <  x<  0.  489.  They 
suggest  that  XQa  is  less  than  0.415  over  the  entire  range  of  temperatures  in 
contrast  to  Fig.  24.  Estimating  X(jQ~  0,40  and  equilibrating  the  partials  of  Zr 
and  N  across  the  a  4  O'  field  yields 

Fn+  =  90,000  -23. 3T  (56) 
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and 


FZr+  +  2RTln2a  =-3400  +  1. 5T 


(57) 


Eq.  (57)  can  be  used  to  calculate  the  vapor  pressure  of  nitrogen  as  a  function 
of  composition  since  (Eqs.  2,  A  1-12) 

T^[x]  -l/2F^  =  l/2RTlnpJ2[x]  =  -FN+  +  RT  lnx(l-2x)-1  (58) 


these  calculations  yields  good  agreement  at  the  low  nitrogen  end  (x  =0.41 5) 
where  the  observed  values  are  one  and  one  half  times  the  calculated  values.  At 
the  high  nitrogen  end  (xw  0,489)  the  observed  value  is  twelve  times  the  calcu¬ 
lated  pressure.  The  results  of  this  calculation  are  given  in  Table  III. 

If  the  high  nitrogen  limit  of  the  ZrN^phase  is  taken  to  be  Xg-j^  -  0.  53 
O' 

as  in  the  TiN  phase  then 

FZr+*  56,200  -  i.25T  (59) 

and 

RTlnCIf  -  -29,800  (60) 

Thus,  the  Zr-N  and  Ti-N  cases  yield  values  of  RTlnfl!  for  the  O' 
phase  intermediate  between  the  carbide  and  oxide  cases.  The  calculations  show 
that  the  free  energies  of  vacancy  formation  decrease  with  increasing  temperature 
for  each  of  the  systems  examined.  Table  IV  contains  a  summary  of  values  ob¬ 
tained  for  the  vacancy  parameter. 
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TABLE  III 


COMPARISON  OF  CALCULATED  VAPOR  PRESSURE  OF 
NITROGEN  OVER  ZrN  (Eq.58)  AND  VALUES 
MEASURED  BY  Klbler,  Lyon,  and  DeSantis*47) 

AS  A  FUNCTION  OF  COMPOSITION  AT  2000°K. 


x  (at.  fr,  N)  1/2  log  (observed)  1/2  log  (Calc.) 

2  2 


0.489 

-2.30 

-2.90 

0.485 

-2.69 

-3.34 

0.480 

-2.97 

-3.67 

0.470 

-3.32 

-3.86 

0.455 

-3.64 

-4.05 

0.440 

-3.92 

-4.15 

0.415 

-4.16 

-4.36 
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TABLE  IV 


CALCULATED  VALUES  OF  THE  VACANCY  PARAMETER  FOR 
SEVERAL  NaCl  TYPE  PHASES 


-RT  In  a/(l-2a) 


PHASE 

(Cal/g.  at 

ZrC 

44, 100 

TiC 

43,850 

TaC 

38,500 

ZrN 

29,800 

TiN 

24, 700 

TIO 

8,200 

Ol  =  number  of  vacant  metal  {or  non-metal) sites  divided  by  the  total  number  of 
sites  at  stoichiometry. 


2a  =  number  of  vacant  metal  (or  non-metal)  sites  divided  by  the  number  of  metal 
(or  non-metal)  sites  at  stoichiometry. 
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3,0  X-Ray  Debye  Temperatures  of  HfC  and  ZrC 

One  of  the  essential  ingredients  in  the  thermodynamic  analysis  of  the 
phase  equilibria  of  a  binary  system  is  a  knowledge  of  the  room  temperature 
entropy  S[  298]  of  a  given  compound  AnBm#  If  one  knows  the  experimental 
value  of  S[  298] ,  one  can  calculate  a  Debye  temperature  parameter  by  one  of 
two  methods: 


S[29«]  = 


6a. 


S[  W8]  »  +  nTi 


m  SD 


f  f  \l/2  A  I 


oA  = 


M 


lxl/2  0 


(61) 

(62) 


B 


The  first  technique  assumes  the  existence  of  some  average  Debye  tempera¬ 
ture  0  which  can  then  be  used  to  calculate  the  specific  heat  or  entropy  by  using 

0 

tabulatod  Debye  functions  Sq[tjt]  .  Electronic  specific  heat  and  Cp-Cv  correc¬ 
tions  are  assumed  negligible  here.  The  second  "two  Debye  0"  approach  is 
based  on  the  empirical  observation  that  two  8's  associated  with  the  individual 
atomic  species  often  yield  a  more  accurate  calculated  specific  heat  curve  when 
compared  with  experimental  data  than  does  the  average  0.  Since  0A  and  9g  are 
not  independent,  either  quantity  is  uniquely  determined  by  S[  298] . 

When  data  is  not  available,  a  direct  measurement  of  (T  or  0^  and  the  use 
of  the  Debye  functions  corresponding  to  either  method  will  yield  specific  heat 
or  entropy  curves.  The  calculated  value  of  S[  298]  compared  with  the  experi¬ 
mental  value  is  considered  to  be  a  relatively  sensitive  test  of  the  accuracy  of 
any  particular  approach.  The  previous  experimental  x-ray  study  (ASD-TR-61- 
445)  of  a  group  of  NaCl  type  compounds  in  which  both  F  and  0^  were  measured 
and  calculated  values  of  S[  298]  compared  with  published  values  was  inconclusive. 
The  small  x-ray  intensity  changes  resulting  from  the  temperature  change  between 
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about  200°C  and  liquid  nitrogen  temperatures  led  to  experimentally  uncertain 
values  of  the  Debye  temperature.  The  present  experimental  program  is  an  at¬ 
tempt  to  extend  the  intensity  measurements  to  a  high  temperature  range  between 
800°C  and  1500°C  relative  to  room  temperature.  The  intensity  changes  for  this 
range  of  temperature  are  quite  large  and  the  resulting  8's  should  be  more 
accurate. 

HfC  and  ZrC  have  been  examined.  Since  good  thermodynamic  data  is  avail¬ 
able  for  ZrC  and  the  compounds  are  assumed  to  be  similar,  the  results  for  ZrC 
can  be  compared  with  the  data  to  check  the  validity  of  the  measurement  while  the 
results  for  HfC  can  be  used  to  predict  specific  heat  and  entropy  curves. 

Thermal  expansion  coefficients  for  the  *wo  carbides  are  also  presented.  In 
the  process  of  obtaining  the  high  temperature  intensity  plots,  the  angular  shift 
of  the  diffraction  peaks  are  also  observed  as  a  function  of  temperature.  Thermal 
expansion  of  the  lattice  constant  is  then  directly  determined. 

The  analysis  of  the  measured  temperature  dependence  of  the  integrated 
x-ray  peak  intensities  from  powder  specimens  of  the  interstitial  compounds  HfC 
and  ZrC  is  based  on  the  following  assumptions: 

1)  The  mean  square  thermal  displacements  of  the  individual  metal  and 
non-metal  atoms  are  effectively  equal  to  a  common  isotropic  dis¬ 
placement. 

2)  The  temperature  dependence  of  the  observed  peaks  can  be  expressed 
in  terms  of  a  Debye  temperature  appropriate  to  an  artificial  mona¬ 
tomic  lattice  of  mass  equal  to  the  average  atomic  mass  of  the  two 
atoms  or  that  of  the  metal  atom. 

The  first  assumption  is  suggested  by  both  experimental  observation  and 
physical  argument.  The  physical  argument  is  based  on  the  fact  that  HfC  and  ZrC 
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are  high  melting  point  compounds.  This  implies  a  relatively  "stiff"  lattice  with 
strong  interatomic  bonding.  Although  different  amplitudes  of  vibration  for  the 


individual  atoms  are  possible  because  a  diatomic  lattice  permits  optical  as  well 
as  acoustical  modes,  the  strong  interatomic  bonding  suggests  that  equivalent 
force  fields  act  upon  each  atom  to  make  the  thermal  amplitudes  approximately 
equal.  Experimentally,  Houska^^  has  measured  the  individual  thermal  ampli¬ 
tudes  of  Zr  and  C  in  ZrC  and  found  the  assumption  of  equal  amplitudes  valid 
within  experimental  uncertainty.  The  expression  for  the  integrated  peak  intensity 


of  a  ZrC  powder  specimen  having  the  NaCl  structure  is: 

4  s.  -BrrCTTT2  -8ff  2  f  t)2 

JJO.T]  -f.  r*i]e  Zr  *f  [^-]e  cL 


where  J  ,  [4» ,  T]  = 


L  [*  .T) 


is  the  normalized  peak  intensity.  Hore 


I±[4>  ,  T]  is  the  observed,  integrated,  peak  intensity  at  Bragg  angle  <1>  at  tem¬ 
perature  T,  p[  hkl]  the  multiplicity  factor  for  peaks  of  Miller  indices  {hkl}  , 

L[<t>]  the  Lor entz  polarization  factor,  and  K  an  instrument  constant.  f„  and 
f  c  are  the  electronic  scattering  factors  for  the  individual  atoms  and  ^ 

and  pc[  T]  2  are  the  corresponding  mean  square  thermal  amplitudes  of  vibra¬ 
tion.  The  +  sign  is  associated  with  all  even  Miller  indices,  the  -  sign  with  all 

J,  +  J_ 

odd  indices;  mixed  indices  are  not  allowed.  If  one  plots  ln(— y, - T.)  and 

J+  "  J-  sin2<t>  Zr 

lnf-jrj - }  against  — ^ -  at  fixed  temperature  then  the  negative  slopes  of  the 

ct  _  _  ^ _ 2  _ 2 

resulting  lines  yield  P^r  and  ^c  respectively.  Because  the  scattering 
power  of  the  metal  atom  is  larger  than  that  of  the  non-metal  atom,  the  determi¬ 
nation  of  pc  2  is  less  certain. 

In  the  present  study,  a  similar  analysis  was  performed  for  HfC  and,  again, 
no  substantial  difference  in  thermal  amplitude  could  be  found. 
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The  second  assumption  that  the  compounds  can  be  treated  as  artificial, 
monatomic,  Debye  solids  can  be  verified  if  the  intensity  data  over  a  large  tem¬ 
perature  range  yields  an  essentially  constant  Debye  temperature  parameter  for 
a  choice  of  either  average  mass  ^  or  metal  atom  mass  MZr  or  M^.  When 
the  metal  mass  is  chosen,  the  model  is  one  of  a  sublattice  of  heavy  atoms  im¬ 
bedded  in  a  carbon  continuum.  The  relationship  between  Debye  0  and  intensity 


I  for  a  specific  peak  (hkl)  is: 

I 


ln(f  )  .  -if— 
o  o 


f[y.  x0] 


(64) 


where  f[y,  xj  «  — { C [ yx  J  -C[xJ  } 


c[x] « 


o  e y  -1 


p  = 


3h2N 


-  a  5.73  x  10"13  c.g. s. 


(65) 

(66) 


2  ,  2  ,  .  2  ,  ,2 

n  =  h  +  k  +1 


y  «= 


e 

t: 


M  , 
M  *  M. , 


e«e 

0  R  0. 


aQ  is  the  lattice  constant,  M  the  atomic  weight,  Nq  Avogadro's  number  and  all 
temperatures  are  expressed  in  °K.  The  function  f[y,  xQ]  has  been  tabulated 
(ASD-TR-61-445  pt.  II)  for  fixed  values  of  y  so  that  high  temperatures  T  are 
preselected  once  the  reference  temperature  Tq  (24°C  in  our  case)  is  given. 
Thermal  expansion  has  the  effect  of  changing  aQ  with  temperature;  however, 
this  correction  was  found  to  be  negligible  over  the  temperature  range  involved. 
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3. 1  Experimental  Results 

The  x-ray  apparatus  consisted  of  a  Picker,  horizontal  diffractometer 
provided  with  a  bent  silicon  crystal  monochromator  and  an  evacuated,  high  tem¬ 
perature  specimen  chamber.  The  carbide  specimens  were  obtained  from  A.  D. 
MacKay,  Inc.  with  reported  purities  of  >99.8a/o.  Both  lattice  parameters  aQ 
were  measured  to  four  figure  accuracy  by  extrapolation  techniques.  For  HfC, 
we  obtain  ao  *  4.  621 A  to  be  compared  with  4.  641 A  reported  for  HfCj  q  and 
4.  61 A  for  the  Hf-HfC  phase  limit^®^,  For  ZrC,  we  obtain  aQ  *  4.695  to  be 
compared  with  *  4.702  for  ZrCg  and  aQ  *4.691  for  ZtCq  55^^*  Both 
materials  lie  on  the  metal  rich  side  of  stoichiometry.  325  mesh  specimen  pow¬ 
der  was  deposited  directly  on  a  tantalum  heating  strip  whose  temperature  was 
monitored  with  a  calibrated  optical  pyrometer  through  a  quartz  window.  A  beryl¬ 
lium  window  provided  access  to  the  specimen  for  monochromated  CuKQ  radiation. 
Conventional  scintillation  counter  and  pulse  height  analyzer  electronics  were  used 
to  display  diffraction  peaks  whose  integrated  intensities  were  measured  graphi¬ 
cally  with  a  planimeter.  Small  specimen  size  and  primary  beam  intensity  loss 
from  crystal  monochromation  led  to  low  intensity  peaks  at  the  higher  Bragg 
angles.  The  improved  signal-to«background  signal  from  monochromation,  how¬ 
ever,  was  sufficient  to  offset  the  difficulty  of  measuring  low  intensity  peaks  and 
relatively  small  intensity  changes  with  temperature. 

The  thermal  expansion  results  are  given  in  TABLE  V.  Three  lines 
each  of  ZrC  and  HfC  were  observed  in  the  temperature  range  800°C  to  1500°C 
at  100°C  intervals.  In  all  cases,  the  angular  shift  with  temperature  was  found 
to  be  linear,  so  that  Of,  the  linear  expansion  coefficient,  could  be  extracted  directly 
from  the  slope  of  each  vs.  T  line.  A  specimen  of  pure  Ta  powder  was  em¬ 
ployed  initially  to  calibrate  the  apparatus  and  the  resulting  data  from  the  Ta 
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powder  indicated  that  extraneous  thermal  displacements  of  the  apparatus  in  this 
temperature  range  would  introduce  errors  not  in  excess  of  5%.  The  expansion 
coefficient, a , is  related  to  the  slope  the  expansion  curve  by: 


a 


TT 

T5D 


cot  q> 


dA<|>  (deg) 

o  3T 


(6?) 


TABLE  V 

EXPANSION  COEFFICIENTS  FOR  HfC  AND  ZrC 

HfC 


line 

d  A<{> 

TT" 

a  (/°C) 

(440) 

70.55 

1.02  x 10"3 

6.29  x  10~6 

(422) 

54.60 

. 540  x 10~3 

6.70  x  10~6 

(420) 

48.19 

. 441 x 10”3 

6.86  x  10"6 

average 

a  «  6. 62  x  10"6/°C 

ZrC 

(440) 

68.11 

1.25  x  10"3 

8.75  x  10“6 

(422) 

53.47 

. 730  x 10”3 

8.45  x  10“6 

(420) 

47. 18 

. 502x  10”3 

8.09  x  10"6 

average  01  =  8.43  x  10-^/°C 

The  expansion  coefficient  for  ZrC,  a  =  8.43  x  10“V°C  can  be  compared 
with  a  measurement  by  Houska^*^  between  754°C  and  1545°C,  namely 
a  a  7.20  x  10‘6/°C. 
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Debye  temperature  results  are  given  in  TABLE  VI.  The  data  was  obtained 
for  HfC  by  measuring  the  temperature  dependence  of  the  (440),  (422)  and  (420) 
lines  under  CuKa  radiation  at  four  temperatures  above  room  temperature, 

24°C.  For  ZrC,  six  temperatures  and  two  lines  (422)  and  (331)  were  used. 

Since  these  compounds  exhibit  cubic  symmetry,  the  9's  obtained  for  any  line 
should  be  the  same.  Consequently  each  number  listed  in  TABLE  VI  is  the  av¬ 
eraged  result  at  that  temperature  of  two  or  three  lines.  The  results  from  line 
to  line  do  not  differ  by  more  than  about  50°K;  however  the  variation  in  0  from 
temperature  to  temperature  is  unexpectedly  large  for  HfC.  Between  830°C  and 
1350°C,both  0  and  decrease  by  about  100°K  whereas  the  0's  for  ZrC  are 
constant  over  this  temperature  range  within  experimental  error.  0  is  expected 
to  decrease  with  increasing  temperature  because  the  atomic  volume  is  increas¬ 
ing.  The  Gruneisen  relation  predicts  a  decrease  in  Debye  temperature  of 
AO  a  SOjdyfTji-Tj)  =  where  y  is  the  Gruneisen  constant.  For  HfC,  Of  * 

7  x  10  ,  9  -  400  and  y  ~  2-3  for  most  metallic  solids.  This  predicts  a  A0  of 

only  about  10°K.  The  temperature  average  of  0  and  0^  are  given  as  440°K 
and  313°K  for  S[  298]  determinations.  Houska/^^  measured  "IF  «  587°K  for  ZrC 
over  the  range  25 °C  to  1680°C,  which  is  considerably  higher  than  our  value  of 
S’  a  435°K,  Because  the  0's  for  ZrC  are  relatively  constant  over  the  temperature 
range  considered,  we  conclude  that  ZrC  behaves  more  like  a  Debye  solid  than 
does  HfC. 

TABLE  VII  contains  the  calculated  room  temperature  entropies  using  both 

IF  in  equation  61  and  9^  in  equation  62.  Since  the  0^,9g  method  gives  better 

agreement  with  the  measured  S[  298]  for  ZrC^®\  we  predict  that  S[  298]  for 

HfC  should  be  closer  to  4.68  than  to  5.92  — -  . 

g.  at.  °K 


-  64  - 


TABLE  VI 


DEBYE  TEMPERATURES  FOR  HfC  AND  ZrC 


Temperature 

F 


®Hf 


HfC 

830°C  907°C 

546°K  439°K 

380°K  323°K 

average  for  F  a  440°K 
average  for  0^  a  313°K, 


1 107°C  1350°C 

382°K  350°K 

282°K  267°K 

Q  a 1205°K 


Temperature 

F 


0 


Zr 


ZrC 

915°C 

1 1 1 5  °C 

1212°C 

470°K 

41 9°K 

427°K 

350°K 

31 5°K 

318°K 

1359°C  1474°C 
430°K  429°K 
326°K  321 °K 


average  for  F  a  435°K 

average  for  a  326°K,  0C  a  897°K 


TABLE  VII 

ROOM  TEMPERATURE  ENTROPIES( - £al  at  298) 

g.  at.  °K 

FOR  HfC  AND  ZrC 


F  method 

ZrC  5.99 

HfC  5.92 


0^,  9g  method 

5.05 

4.68 


measured 


(50) 


4.01 
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4.  0  Estimation  of  the  Free  Energy  of  Formation  of  Hafnium  Carbide 

(J 

The  free  energy  of  formation  of  stoichiometric  hafnium  carbide,  AF  [  T] , 
in  units  of  cal/g.at,  can  be  approximated  as  follows: 

AF^tT]  =  AHff[0°K]  +  AF^[T]  (68) 

InEq.  (68),  AH^[0°K]  is  the  enthalpy  of  formation  at  0°K,  while  AF^[T] 
is  defined  by  Eq,  (69), 

AF^[T]  =  j^ACpdT  -  T  j^T"1  ACpdT  (69) 

where  ACp  is  the  specific  heat  of  the  stoichiometric  compound  in  cal/g.  at.  °K 
minus  one  half  the  sum  of  specific  heats  of  hafnium  metal  and  graphite  in 
cal/g.  at.  °K. 

(J 

A  preliminary  value  for  AH  [  298°K]  has  been  obtained  by  K.  K.  Kelly  by 
using  combustion  calorimetry.  This  value  (quoted  on  p.  63  of  15  June  1963  report 
on  contract  AF33(657)-8223,  H.  L.  Schick  et  al  AVCO  Corp,  ,  Wilmington,  Mass.), 
which  is  subject  to  modification  by  Kelley,  has  been  reported  as  -29.  5  kcal/g.  at. 
or  -59.0  kcal/mol.  Since  AH^  will  not  vary  by  more  than  200  cal/mol  between 
0°K  and  298°K,  the  latter  value  can  be  safely  adapted  for  AH*\  0°K]  i.  e. 


AH^[0°]  *  -29.5  kcal/g.  at. 

(70) 

2AH°[0°K]  £  -59.0  kcal/mol 

(71) 

Consequently,  specification  of  AF^[T]  requires  that  AF^f[T]  be  evaluated. 
Since  no  low  temperature  specific  heat  data  are  available  for  HfC,  the 
x-ray  results  obtained  in  section  j.O  of  this  report  will  be  used  to  estimate 

Crr  , 

AFp[T]  .  This  procedure  requires  considering  the  vibrational,  electronic  spe¬ 
cific  heat,  and  phase  change  components  of  hafnium,  carbon,  and  hafnium  car¬ 
bide  to  AF^[  T]  . 
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4.  1  Hafnium 


Walcott^ ^  has  measured  the  specific  heat  of  hafnium  between  4°K 
and  300°K.  He  suggests,  0^  =  200°K,  and  y^  =  6.3  x  10~4  cal/g.at.  °K2, 
Hafnium  undergoes  a  h.c.p.  -*b.c.c.  transition  on  heating  above  2033°K.  We 
can  approximate  AS  ,  i.e.,  the  entropy  change  for  the  h.c,p,-*b,c.c.  tran¬ 
sition  as  0.9  cal/g.at.°K  on  the  basis  of  earlier  work  on  Ti  and  Zr^.  The 
melting  point  of  Hf  is  2495°K,  and  we  approximate  A L  =  2. 1  cal/g.at. °K. 

4. 2  Graphite 

The  free  energy  function  for  the  graphitic  form  of  carbon  1b  tabulated 
by  Janaf  and  by  Elliot  and  Gleiser^2^.  These  data  can  be  used  to  compute  F^  [T] . 
4.  3  Hafnium  Carbide 

Westrum^50^  has  measured  the  specific  heat  of  ZrC  between  5°K 

and  350°K.  Plotting  his  low  temperature  data  (5-25°K)  as  Cp/T  vs  ’f2  yields 

a  very  low  value  of  y  for  ZrC  (i.e.,  y  Z  0.5  x  10  4  cal/g.at.°K2),  Similar 

values  have  been  noted  for  the  refractory  borides  ZrB2,  TiB2  and  HfB2>  On 

the  basis  of  section  3.0,  we  approximate  0^  =  313°K  and  0^  =  1205°K  for  HfC. 

0 .  . 

On  the  basis  of  this  information  AFT  L T J  can  be  computed  as 
follows:  For  hafnium,  we  evaluate  l T]  by  using  Eq.  (72) 

F°f  IT]  =  Ft-^-]  +  1O‘4(3R02P-T  Ut-4-J)  -  ^Hf1"2/2  +  [t]  4Af|^1t] 

where 

P  =  T2/2  £'2  -  (l/ 40)  In  (1  +  20  T2/0  2  ) 

0  '0 

and  F[  ]  and  U[—  ]  are  the  Debye  free  energy  and  energy  functions.  These 
equations  derive  from  the  approximation 

Cp  =  Cy  [^r]  (14  10"4T)  4  yT 

$ -  o(52) 

Recent  high  temperature  measurements  (1200-1000  K)'  on  NbC  and  TaC 

offer  further  support  for  this  approximation. 


where  6  =  200°K,  and  y  =  6.3  x  10  4  cal/g.at  °K2.  The  first  term  on 

the  right  of  Eq,  (72)  is  the  vibrational  free  energy  term^  (at  constant 

volume),  the  second  term  arises  from  the  (C  -  C  )  correction  which  is 

P  ^ 

A  A 

approximated  by  10"’T  Cy  [  -y-]  in  Eq.  (74).  The  third  term  is  the 
electronic  specific  heat  contribution  to  the  free  energy.  The  fourth  term, 
which  is  added  only  for  T>2033°K  is  given  by 

AFH£“^  S  -  0.9  (T-2033°K)  cal/g.at.  (75) 

while  the  fifth  term  which  is  added  only  for  T>2495°K  is 

ArJ*  S  -2.1  (T-2495°K)  cal/g.at.  (76) 

Or 

A  similar  procedure  can  be  used  for  to  compute  F^,[  T  ]  for  HfC. 
by  using  an  expression  similar  to  Eq,  (72)  with  0^  =  1205  °K  and  y  _  0 
i.  e. 

ZFT°[T]  »  F[-^]  +  +  3xl0'4R((313)2  P[313]  +  (1205)2  P 1 1205] ) 

-  10’4T(U[^y]  +  U[-^]  )  cal/mol  (77) 

Since 

2AFxg  It]  =  2FFg[T]  -  F^f[T]  -  f'J.  [  T  ]  (78) 

Performing  the  numerical  calculations  of  Eqs.  (72)  to  (78)  and 
combining  with  Eqs.  (71)  and  (68)  yields  the  values  for  the  free  energy  of 
formation  which  are  tabulated  in  Table  VIII.  These  values  have  been  used 
in  connection  with  the  calculations  of  vapor  pressure  and  phase  equilibria 
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presented  in  section  2.0,  (See  pp.  22,  28,  29,  36,  and  37.)  In  view  of  the 
agreement  shown  in  Figs.  18  and  19,  it  appears  that  the  values  in  Table  VIII 
are  accurate  to  within  2  kcal/mol  or  better. 


TABLE  VIII 

ESTIMATED  VALUES  FOR  THE  FREE  ENERGY  OF  FORMATION 
OF  STOICHIOMETRIC  HAFNIUM  CARBIDE 
(2AF°[0.5,  T]  in  kcal/mol) 

2AH<J[0.5,  0°K]  2  2AH°  [0.5,  298°K]  s  -59. 0  kcal/mol 


T°K 

2AFff[T]  (kcal/mol) 

1400 

-55.8 

1600 

-55.4 

1800 

-55.0 

2000 

-54.6 

2200 

-54.2 

2400 

-53.6 

2600 

-53.  1 

2800 

-52. 1 

3000 

-51.2 

3200 

-50.2 

3400 

-49.3 

3600 

-48.3 

3800 

-47.3 

4000 

-46.3 
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APPENDIX  A1 


EXPLICIT  FORMULATION  OF  THE  COMPOSITIONAL 
DEPENDENCE  OF  THE  FREE  ENERGY  OF 
NON-STOICHIOMETRIC  PHASES 

In  considering  the  thermodynamic  properties  of  refractory  oxide,  car¬ 
bide,  and  nitride  compounds  it  is  quite  apparent  that  many  of  these  compounds 
are  not  line  compounds.  In  reality  they  are  found  to  exist  over  a  range  of  com¬ 
positions  at  elevated  temperatures  and  even  at  room  temperature.  Recently 
these  deviations  from  stoichiometry  have  been  recognized  as  contributing  fac¬ 
tors  in  controlling  many  of  the  important  properties  of  these  compounds.  In 
order  to  proceod  with  the  analysis  of  the  thermodynamic  properties  of  metal- 
nonmetal  systems  it  is  necessary  to  derive  equations  which  can  be  used  to 
explicitly  define  the  thermodynamic  properties  of  these  nonstoichiometric  com¬ 
pounds.  This  conclusion  is  apparent  if  one  considers  TiO,  ZrC,  TiN,  or  TaC 
as  examples. 

Schottky  and  Wagner^^  have  formulated  a  generalized  thermodynamic 

model  of  non-stoichiometric  compounds.  A  detailed  description  of  this  general 

(191 

model  is  given  by  Wagner'  '  in  which  one  gram  atom  of  a  compound  A.,  .B 

v  t  —  xj  x 

is  considered.  In  the  present  notation  x  is  the  atomic  fraction  of  element  B 
(rather  than  Wagner’s  symbol  x^).  The  general  model  allows  for  two  sublattice, 
an  A  sublattice  and  a  B  sublattice,  which  contain  both  A  and  B  vacancies,  A 
and  B  interstitials,  A  atoms  on  B  lattice  sites  and  B  atoms  on  A  lattice  sites. 
Restricting  consideration  to  the  interstitial  type  compounds  of  interest  i.  e. , 
where  the  B  atom  is  much  smaller  than  the  A  atom  and  the  B  sublattice  fits 
within  the  interstices  of  the  A  sublattice,  the  general  model  can  be  restricted. 
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As  a  first  approximation  we  consider  the  case  where  the  only  defects  present 
are  vacant  A  sites  and  vacant  B  sites. 

We  consider  a  compound  A^_xjBx  which  has  the  6  crystal  structure  with 

a  composition  range  about  x  .  Thus  the  compound  is  ideally  An  .B  ,  but 

0  '1“xo'  xo 

can  exist  for  values  of  x  greater  than  xq  (corresponding  to  an  excess  of  va¬ 
cant  A  sites)  or  less  than  xq  (corresponding  to  and  excess  of  vacant  B  sites). 
Since  we  consider  a  gram  atom  of  compound  then  the  number  of  A  atoms  is 
given  by  Eq.  (Al-I)  as 

N(l-x)  »  (l-xo)Ng-NA+  (Al-I) 

where  N  is  Avogadro  ’s  number,  Ns  is  the  total  number  of  lattice  sites  (A  sites 
plus  B  sites)  and  is  the  number  of  A  lattice  sites  which  are  vacant.  The 
number  of  B  atoms  is  given  by  Eq.  (Al-2)  as 

Nx  »  xQ  Ng  -  Nb+  (Al  -2) 

where  is  the  number  of  vacant  B  lattice  sites.  The  ratio  of  the  total  num¬ 

ber  of  sites  Ng  to  the  total  number  of  atoms,  Ns/N,  is  set  equal  to  a  param¬ 
eter  y.  Hence  the  fraction  of  vacant  sites  is 

(N8-N)/Ng  «(l-l/y)  (Al -3) 

The  free  energy  of  the  6  phase  as  a  function  of  temperature  and  composi¬ 
tion  (reference  (19)  equations  (3-1) -(3-8))  is  given  as 

F6[T,x]  =  ( 1  -x)  [  T ]  +  xF°[T]  +  (NA+/N)FA+[  T] 

+(N8/N)(F*  [Tl  -(l-xo)F°  [T]  -x0F°[T])+(NB+/N)FB+[T] 

+  kT  (NA+ln(NA+/(l-xo)Ns)  +  *B+  ln(NB+/xoNs) 
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+(  (l-x0)N8-NA+)ln((l“Xo> 


+  (xoN8-NB+)ln(xoN8-NB+)/xoN8>  {A1-4) 

These  equations  apply  for  the  case  where  f^[0oK]  and  E^)[0°K],  which  are 
the  free  energies  of  one  gram  atom  of  A  and  B  respectively  at  one  atmosphere 
pressure,  are  equal  to  zero  at  0°K,  and  all  of  the  free  energies  are  understood 
to  be  temperature  dependent.  In  these  equations  FA+  and  Fg+  are  the  free 
energies  of  a  gram  atom  of  A  and  B  vacancies. 

For  a  given  temperature,  pressure  and  composition  i. e.  (N^  and  Ng  are 
constant),  the  free  energy  F  is  a  minimum  whent 

(Al-5) 

-RTlna  =  S*6[xo]  -  (1-*0>(FA  “  ^A+  +  RT In (1  -xo))-xo(Ej^» Ffi+  +  RT lnxQ) 
and 

lnff  s  (l-xQ)  ln{(l-xQ)  y  +  x-l)/y  +  xq  In  (xQy-x)/y  (Al-6) 

Substitution  of  these  expressions  into  (Al-4)  yields  F ^  T ,  x]  given  in 
TABLE  Al-I.  The  partial  molar  free  energies  ^  and  Tg  are  then  derived  by 
applying  the  general  relations, 

FA“r--S  "O  If 

The  equations  given  in  Table  Al-I  should  apply  for  situations  where  x  is  not 
very  different  from  xq,  and  the  fraction  of  vacant  sites,  (Ng-N)/Ng,  is  small 
for  x=xo  i,  e.,  (y— 1)  is  small  at  xsxq.  Under  these  conditions, 

^  F*  [x  ]  is  the  free  energy  of  one  gram  atom  of  vacancy  free  stoichiometric- 
compound 
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TABLE  Al-I 


THERMODYNAMIC  EQUATIONS  FOR  ONE  CRAM -ATOM  OF 

NONSTOICH IOME T RIC  COMPOUND 

A,.  ,  B  WHICH  HAS  THE  6  STRUCTURE 

(1-x  )  x 
'  o'  o 

x  =  atom  fraction  B  N^+  =  number  of  vacant.  A  sites 

N  s:  Avogadro’s  Number  NB+  =  number  of  vacant  B  sites 

N  »  total  number  of  lattice  sites  y  =  N_/N,  y>  1 

General  Case 

N(l-x)  =(l-xo)  Na- Na+  Nx  =  xqNs  -  NB+ 

Ff  a{l-x)  F°  +  xF°  +  y(Pl<[x0]  -(l-xo)F^-xoF“  )  1  ((l-xo)y-(l-x))FA+ 
+{xQy-x)FB+  +  RT((y(l-xo)  +  x-l)  ln((l-xQ)y+ x-l)/y{l-XQ) 

+(xQy-x)  ln(xoy-x)/xoy  +  xlnx/xQy 
+  <l-x)  ln(l-x)/(l-xo)y) 

F|  =  F°  -  Fb+  +  RT  lnx/(xoy-x) 

FA  =  FA  “  FA+  +  RTln{l-x)/{(l-xo;.  y-h  x-l) 
where 

In  a  =  ( 1  -xQ)  ln((l  -xQ)y  +  x-i)/ y  +  xq  in  {xQy-x)/y 

and 

-RT  In  a  =f|[xo]  -(l-xo){F°-FA+  +  RTln(l-xo))  -xq(F°-Fb+  +  RTlnxQ) 
hence 

,  ,  -x  (x  -1) 

AF°  =  A  F*  +  RT  In  ( 1  -  a  xq  °(1-xo)  °  ) 
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TABLE  Al-I  (Continued) 


when  x  =x 

o 

X  ]  -X  X  1  “X 

a-xo°(l-xo)  °{l-l/y)axo  °(l-x)  °(Ns-N)/Na 

and 

NA+°xo<VN>5 

Special  Case 
xQ  =  0.  5  ,  Aq>  5Bq>  5 

N{l-x)  r  0.5  Ns-Na+  ,  Nx=0.5N8-Nb+ 

F6a(l-x)F^+  xF°  +  y(Ft?  [  0.  5]  -  0.  5  F°  -  0.  5  FB°)  +  (0.  5  y-x)  FR+ 

+(0.5y+x-l)FA+  +  RT({0.5y+x-l)  ln(0.  5y+  x-l)/0.  5y+  xlnx/0.  5y 
+(0.  5 y-x)  ln(0,  5y-x)/0,  5y+(l-x)  ln(l-x)/0.  5y) 

^A  "  FA-FA+  +  RTln{l-x)/(0.5y-(l-x)) 

Fb  =  F°  -  F„+  +  RT  lnx/(0,  5y-x) 

where  -RT  lnCH  =  Ffc6  [  0.  5}  -  0.  5(FA  +  F°  -  FA+ -  FB+)  -  RT  InO.  5 
and  y  =  (l-4a2)"1(l  +(l-4x(l-x)(l-4a2))1'/2) 

when  x=  0.5,  y*l/(i-2a),  and  N../N  »Nn,/N  s«aO,5(N  -N)/N 

A  T  8  XJ  T  9  S  53 

AF6[0.5]  =AF#(0.5]  +  RTln(l-2a) 
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F>:<  t  x0^  "  ^ 1  "xo)  FA  "  xoFB 

is  the  free  energy  of  formation  of  one  gram  atom  of  the  stoichiometric  com¬ 
pound.  For  a  fixed  number  of  A  and  B  atoms  (i.  e.  a  fixed  value  of  x),  F^ 
given  in  Table  Al-I  is  a  minimum  when  y  takes  on  a  specific  value.  This  value 
of  y  is  specified  by  means  of  the  parameter  a  in  terms  of  F^+ ,  Fg+,  and  the 
free  energy  of  formation  of  the  ideally  ordered  compound  [xq]  »  (for  which 

y  =  !)• 

The  equations  for  the  case  of  a  general  value  of  xq  are  simplified  consider¬ 
ably  if  we  consider  a  special  case  when  xa0.5  i.  e.,  TIO,  TiN,  ZrC  etc. 

In  this  case  we  see  that  a  is  equal  to  the  fraction  of  vacant  A  lattice  sites 
(and  the  fraction  of  vacant  B  lattice  sites)  at  x=0,5.  The  equations  for  the  par¬ 
tial  molar  free  energies  of  A  and  B  are  complicated  by  the  fact  that  y  is  a 
complex  function  of  x  even  for  the  special  case  of  xq  =0.  5.  However,  if  we 
recognize  that  a  (which  is  the  fraction  of  vacant  A  or  B  sites  at  x*D,5)  is 
small  then  the  expression  for  y  can  be  simplified.  The  usual  values  for  a  are 
considerably  less  than  1%.  If  we  restrict  consideration  to  the  situation  where 
or  is  very  much  less  than  unity  then 


(N  /N)  *ysl*[(l  -2x)  (Al-8) 

8 

since  y  must  be  greater  than  unity  the  plus  sign  refers  to  x<0.  5  and  the  minus 
sign  to  x>0.5  in  eq.  Al-8,  Under  these  conditions  the  following  expressions 
result:  at  x  =  0.  5 


FA  =  FA  -  fa  +  RT  ln  d  -Za  >/2Q; 


(A  1-9) 


=  FB  -  FB  +  RTln(l-2a)/2<* 


(Al-10) 
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Eqs,  (Al-9  andAl-10)  apply  even  when  a  is  not  small. 

2  2 

When  x<  0.5  (and  a  /(l-2x)  is  very  small)  then 

=  FA_IA+  +  RTln(l-2x)/4(l-x)a2  (Al-U) 

=  Fg  -  Fg+  +  RT  ).nx/(l  *2x)  (Al-12) 

while  for  x  >  0,  5 

*A  =  FA  ‘  FA+  +  RTln(l-x)/(Zx-l)  (Al-13) 

Tg  =  Fg  -  Fb+  +  RTln(2x-l)/4xa2  (Al-14) 

If  aA[x]  and  aA[x0]  are  the  activities  of  A  at  a  composition  x  and  at 
xQ  a  0. 5  then  for  x<0.  5  the  ratio  of  activities  (or  vapor  pressures)  is 

aA[xl/aA[xQ]  =  (l-2x)/2(l-x)or  (Al-15) 

and 

ag(  X1  / ag[  XQ]  =  2  O  x/(l  -2x)  (A  1-1 6) 


Thus  the  activity  of  A  increases  while  the  activity  of  B  decreases  with 
increasing  deviation  from  stoichiometry.  On  the  other  hand,  when  x>xq. 

aAt  x] /aA^  Xo]  b2«(1-x)/(2x-D  (Al-17) 

and 

aB[x]/aB[xo)  =  (2x-l)/2  xor  (Al-1^) 

so  that  the  activity  of  A  decreases  and  the  activity  of  B  increases  with  in¬ 
creased  deviation  from  stoichiometry.  This  is  the  general  behavior  to  be  ex¬ 
pected  for  such  compounds. 

At  this  point  it  is  worthwhile  introducing  a  simplification  by  inserting  the 
relation  between  a,  F^+,  and  FRJ_  into  the  equations  for  and  If  we 


B  + 
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let  AFj|  x^.T]  be  the  free  energy  of  formation  of  one  gram  atom  of  the  com¬ 
pound  i,  e.  * 

AFfUo,T]  .  F5[xJ  -(1-xo)f£~xoFb°  (Al-19) 

then 

-2  RT  In  a  h  2.  A  Ff  [  xo>  X]  +  F^+  +  FB+  +  2  RT  In  2  (A1  -20) 

Hence  from  eqs,  (Al-11-15)  we  find,  (with  xq  *«  1/2) 
for  x<  0.  5 

VA  sFA+  2aF£[1/2,T]  +  Fb+  +  RT  In  (1  -2x)/(l»x)  (Al-21) 

“  %°"  FB+  +  RT  lnx/(l"2x)  (Al-22) 

and  c 

^A  s  FA  "  FA+  +  RTln(l-x)/(2x-l) 

FB  a  F°  +  2AFf[  1/2,  T]  +  Fa+  +  RT  In  (2x»l)/x  (Al-23) 

for  x>  0.  5. 

These  expressions  can  be  used  to  calculate  the  compositional  variation 
of  the  enthalpy  of  formation  of  the  6  phase  at  298°K.  If  O'  is  small  (which  is 
probably  likely  at  low  temperatures)  then  the  enthalpy  of  formation,  AH£,  is 
given  by 

AHf  =  H6-(1-x)  H^-xHg0  =  yAH:;:  +  (y/2  +  x-l)HA+  +  (y/2-x)HB+  (Al-24) 

where  AH*  is  the  enthalpy  of  formation  at  stoichiometry.  For  x  less  than  0.  5 

Note  that  these  approximations  apply  only  when  F^[xq]  =  F^[xq]  or  when  O' 
is  very  small.  If  this  is  not  so,  then  a  more  exact  expression,  namely 
F^xJ  =  fI[xq]  +  RT  In ( 1  -2a)  results. 
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(i.  e.  y-2(l-x))  then 

AHf -AH*  *  (l-2x)(AH*  +  Hb+) 
while  when  x  is  greater  than  0.5(y-2x) 

AHj-AH#*  (2x»l)(AH^t  +  HA+) 


(Al-25) 


(Al-26) 


Application  of  these  equations  for  the  case  where  xq  is  different  from 
l/2  in  the  exact  form  is  complicated  by  the  fact  that  the  solution  for  y  in  terms 
of  a  involves  cubic  or  higher  order  equations,  However  some  simplifications  can 
be  applied  if  a  is  very  small. 

Recasting  Eq.  Al-6  as: 

(ay)^^  a  (xoy-x){l -xQ)y  +  x »  l)“^“*/xo^  (Al-2?) 

indicates  that  when  a  is  small,  y  -  x/xq  for  x>xq  and  y*  (1  »x)/(l  -xq)  for 

x<xq.  These  approximations  can  be  used  to  write  the  analogues  of  Eqs.  Al-9- 

14  for  x  4  l/2:  thus  at  x  =  x  (for  all  values  of  O') 
o  o 


rA 

when  a  is  small  andx<x. 


F°-Fa+  +  RT  In  (a  "1xo°  (l-xo)'l“xo)  -1) 

(Al-28) 

FB°"FB+  +  RTln(a''1x*°(l-x0){1“Xo)-l) 

(Al  -29) 

d-xj 


^f*FA"FA+  +  RTln{l-x)  °(xo-x)  °a_1 

TB  “FB-FB4  +  RTlnxd-x^fx^xr1 


and  when  x>x 


*  Fa°-Fa+  +  RTlnxJl-x)(x-xo) 


-1 


LA  A+ 


(Al-30) 

(A  1—31) 

(Al-32) 
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(Al-33) 


s.  \  (x  -1)  x  (1-x  )  , 

*  FB°"  FB+  +  xo  RTlnx  xo  (x“xo)  ° 


Similarly  the  compositional  dependence  of  the  enthalpy  of  formation  is 
given  by  equations  (A1 -34-35)  as: 


“1 1 


AHf~AH#  «  (xo-x)(l-xo)  ‘(AH,  +  Hb+) 


(Al-34) 


for  x<  xq  ,  and 

AHf-AH*  »  xo"1(x-xo)(AH«<  +  hA+)  (Al-35) 

for  x>  x 

o 
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APPENDIX  A2 


APPROXIMATION  OF  THE  EXCESS  FREE  ENERGY  OF 
MIXING  OF  INTERSTITIAL  SOLID  SOLUTIONS 

(4) 

Kaufman,  Radcliffe,  and  Cohen'  '  have  formulated  the  compositional  and 
temperature  dependences  of  an  a  interstitial  solid  solution  of  element  C  in 
element  A  in  which  there  exist  n  interstitial  C  sites  per  A  site.  In  this  model 
z-1  of  the  interstitial  sites  surrounding  a  given  "G-site"  are  excluded  from  oc¬ 
cupancy  (by  other  C  atoms)  if  it  is  filled  by  a  C  atom.  In  this  situation,  the  a 
solution  terminates  at  a  composition  corresponding  to  x»n(n+z)  *  corresponding 
to  a  compound  ACn^z  .  The  free  energy  of  a  gram  atom  of  such  a  solution  is  de¬ 
fined  by  Eqs.  5-7  of  reference  (4).  Moreover,  the  partial  molar  free  energies 
are  given  in  terms  of  Eqs.  (7-8)  of  reference  (4),  In  order  to  complete  the  ex¬ 
plicit  formulation  of  the  free  energy  of  the  Of  solid  solution,  it  is  necessary  to 
define  the  excess  free  energy  of  mixing  Fg  and  the  free  energy  of  the  terminal 
compound  AC^j^  . 

Asa  first  approximation,  0°K  enthalpy  contributions  and  changes  in  vibra¬ 
tional  free  energy  are  considered  as  the  primary  contributions  to  these  terms. 
Restricting  consideration  to  high  temperatures,  the  vibrational  free  energy  can 
be  approximated  by  3RTln6/T  where  9  is  the  Debye  temperature.  On  the 

basis  of  these  assumptions  the  free  energy  of  a  gram  atom  of  AC  /  is: 

n/  z 


0, 


*a 


(A2-1) 


FAC  ,  =HAC  /  [°°K1  +  3RT  {(l.-£-)ln-^-  + 


yn/z 


J  n/ t. 


In  Eq.  (A2-1)  9^  and  6^,°  are  the  individual  Debye  Temperatures  of  A 
atoms  and  C  atoms  in  the  compound  d ACn^7  .  According  to  the  Two-Debye- 
0  method^ 
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0c*6a“a2“c/2 


thus 


0. 


*®  -Cma2mc1/2 


(A2-2) 


(A2-3) 


In  Eq.  (A2-D.H? 


. .  _  [  0°K]  is  the  heat  of  formation  of  one  gram  atom  of 

n/  z 

aACn^z  from  (1-^j)  gram  atoms  of  a  A  and  gram  atoms  of  the  SQ  form 
of  C  at  0°K.  Since  the  enthalpy  of  a  A  and  the  SQ  form  of  C  at  0°K  are  both 
zero  at  0°K,  then  in  the  spirit  of  the  present  approximation 

Fj[T]  s  3  RT  In  0^/T 

cpa 


(A2-4) 


at  high  temperatures,  where  0^  is  the  Debye  temperature  of  pure  aA. 

Definition  of  requires  a  division  of  Fg  into  two  terms,  a  0°K  teri 
and  a  temperature  dependent  term.  The  simplest  form  of  the  0°K  term  is 
Ax(l-x(n+z)n”1)  which  is  zero  at  x«0  and  xsn(n+?.)”*.  If  the  Debye  temperature 
of  A  and  C  are  related  by  Eq.  (A2-2)  then; 


n  J 


F®[x,T]  =  Ax  [  1  -  -—x^  3  RT  {  (1  -x)  In 


+  xlnH*a — j  /  ^a) 

'  C  /  v  A  / 


(A2-5) 

z/n) 


a  ot 

where  0^  [  x]  and  0£,[x]  are  the  Debye  temperatures  of  elements  A  and  C  in 
the  a  solution  at  composition  x,  and  A  is  a  constant.  On  the  basis  of  this  defi¬ 
nition,  the  chemical  potentials  of  elements  A  and  C,  and  ,  in  the  a  solid 
solution  can  be  defined  as  (Eqs.  5-9  of  ref.  (4))  follows: 

(A2-6) 


I?  =  F?  +  In  -n--*Ln+?>  +  Ax2  +  3  RT  ln  -£ 


A 

and 


"n(l-x) 


oa 


A1 

~a t 


no a 


to'jqtSRTln-f-tRTtos^sSx 

n/z  ' 


(A2-7) 


+A  (  1 -2x 


tSti  t  l2i£!)  ,  3  Rl  /ln  t  S!-l 

n  n  j  \  qoQ  Qa  r dx 


e"[x] 
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Consequently,  this  formulation  specifies  the  free  energy  of  the  interstitial 

solution  in  terms  of  n,  which  can  be  evaluated  from  a  knowledge  of  the  crystal 

structure  (i.  e.  n  =  1  for  h.  c.p.  and  f,  c.  c.  lattices  and  3  for  b.  c,  c,  lattices); 

the  compositional  dependence  of  the  Debye  temperatures  A  and  C,  which  can 

(5) 

be  estimated  from  a  knowledge  of  the  melting  point  and  volume  ,  the  exclu¬ 
sion  parameter  z,  which  is  equal  to  unity  for  an  Ideal  interstitial  solution; 

and  the  parameters  A  and  H^_  [0°K], 

n/z 

The  enthalpy  of  formation  of  one  gram  atom  of  solution  from  (1-x)  gram 
atoms  of  Aa  and  x  gram  atoms  of  the  SQ  form  of  C(the  stable  form  at  0°K  and 
one  atmosphere  pressure)  is, 

£W,[xf  0°K]  ■  n^xfn+z)  H^r  [0°K]  +  Ax(l -n"*x(n+z)) 


n/z 


(A2-8) 
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TABLE  OF  SYMBOLS 


Several  symbols  used  in  this  paper  have  more  than  one  connotation.  How¬ 
ever  in  the  few  cases  where  two  meanings  are  possible  the  proper  choice  is 
apparent  from  the  text.  It  should  be  noted  that  brackets  are  used  exclusively 
to  denote  functional  relations.  Thus  V^x,  T,P]  is  the  volume  of  the  O'  phase 
which  is  a  function  of  composition,  x,  temperature  and  pressure.  As  is  indi¬ 
cated  in  the  text,  all  vapor  pressure  and  free  energy  terms  are  temperature 
dependent  and  free  energies  are  discussed  in  units  of  calories  per  gram  atom. 


B 

C 

u 

P 


7 

L 


Me 

T7 

Me 

T* 


X 


xj3d 

x/3L 

xCT/3 


a  pure  element  B. 
a  pure  element  C. 

a  phase  having  the  sodium  chloride  crystal  structure, 
a  phase  having  the  b.  c.  c.  crystal  structure, 
cie  "graphite"  crystal  structure, 
the  liquid  phase. 

the  melting  point  of  the  form  of  a  metal, 
the  melting  point  of  the  y  form  of  a  metal, 
the  melting  point  of  the  O'  phase. 

atomic  fraction  of  non-metal  (i.  e,  carbon,  oxygen,  nitrogen), 

the  eutectic  temperature  corresponding  to  a  eutectic  composition  which 
is  less  than  x  » 0.  5. 

the  composition  of  the  0/0+ O’  phase  boundary, 
the  composition  of  the  0/0+ L  phase  boundary, 
the  composition  of  the  d/(3+  0  phase  boundary. 
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xLff  the  comPoa*tion  °f  the  L/L+  d  phase  boundary  for  x^g-  less  than  0.  5, 

the  composition  of  the  d/L.  +  d  phase  boundary  for  Xg^  less  than  0.5. 

xLd  the  cornPos*tion  of  the  L/L+  d  phase  boundary  for  x.  g  greater  than 

*v  or  equal  to  0.5. 


xdL 

if 

A 

F°[x] 


K  x 

A  xa/3 


the  composition  of  the  d/L+  d  phase  boundary  for  3t-,  greater  than 
or  equal  to  0.  5. 

the  free  energy  of  the  (3  form  of  A  which  is  (understood  to  be)  a  func¬ 
tion  of  temperature. 

the  integral  free  energy  of  the  d  phase  which  is  a  function  of  compo¬ 
sition  (and  temperature). 

the  partial  molar  free  energy  of  A  in  the  d  phase  which  is  a  function 
of  composition  and  temperature. 

the  partial  molar  free  energy  of  A  in  the  d  phase  evaluated  at  a  com¬ 
position  corresponding  to  the  boundary  of  the  d/d+  (3  phase  field. 


the  free  energy  of  the  graphite,  y ,  form  of  carbon, 
the  vapor  pressure  of  the  |3  form  of  pure  Zr, 

the  vapor  pressure  of  Zr  over  the  d  phase  which  is  a  function  of  com¬ 
position  (and  temperature). 

the  free  energy  of  formation  of  the  d  phase  at  stoichiometry  (cal/g.at. ), 


^Me+  the  *reQ  ®nergy  of  formation  of  Me  vacancies. 

Fr^  the  free  energy  of  formation  of  C  vacancies. 


a  the  phase  Of,  or  the  vacancy  parameter  which  is  proportional  to  the 

fraction  of  vacant  lattice  sites  at  stoichiometry  (See  Table  AI  for 
mathematical  definition). 

AFp  the  difference  in  free  energy  between  the  6  (diamond  cubic)  and  y 

(graphitic)  forms  of  carbon. 


•y . 

AF '  the  difference  in  free  energy  between  the  y  and  liquid  form  of  carbon. 

s/ 

-y,-* 

AS^i  the  difference  in  entropy  between  the  y  and  liquid  form  of  carbon. 


N 


Avogadro's  number. 


k  Boltzmann^  constant. 
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R 


AF, 


L-0 


Me 


AH  [  0UK] 


G 

XC 

1/2FoG2 

H^.i0[  0°K] 


AF, 


a-'(i 


Ti 


,00! 

Zr 


N 


N, 


A+ 


F°[T,x] 


aA[x] 


a  intersti¬ 
tial  solid 
solution 

n 


z  - 1 


-01 


the  gas  constant  =  1.9873  cal/g.at.  =  Nk. 

the  difference  in  free  energy  between  the  /3  and  liquid  forms  of  a  metal. 

the  enthalpy  of  formation  of  one  gram  atom  of  CT  phase  of  stoichiometric 
composition  at  0°K. 

the  rate  of  Langmuir  evaporation  (mols/cm^sec. ) . 

the  composition  of  congruent  vaporization  (i.  e.  =x^[  T] ). 

the  free  energy  of  one  half  mol  (i.  e,  one  gram  atom),  of  oxygen  gas 
at  atmospheric  pressure. 

the  Debye  temperature  of  Ti  in  the  0!  solid  solution  (which  is  a  func¬ 
tion  of  composition). 

the  enthalpy  of  the  Cl  form  of  TiO  at  0°K. 

element  A,  or  the  enthalpy  of  mixing  parameter  for  a  interstitial 
solid  solutions  (see  Eq.  A2-5). 

the  free  energy  difference  between  the  a  and  /3  forms  of  titanium. 

the  Debye  temperature  of  pure  a  Zr  i.  e,  0£r[x]  at  x"0, 

the  total  number  of  lattice  sites, 

the  number  of  A  lattice  sites  which  are  vacant. 

the  ratio  of  the  total  number  of  sites  to  the  total  number  of  atoms, 

N/N. 

the  f~ee  energy  of  the  6  phase  which  is  a  function  of  temperature 
and  pressure. 

the  free  energy  of  one  gram  atom  of  vacancy  free  stoichiometric 
compound. 

the  activity  of  element  which  is  a  function  of  composition  and  under¬ 
stood  to  be  a  function  of  temperature. 

the  interstitial  solid  solution  of  element  C  in  element  A. 

the  number  of  interstitial  C  sites  per  A  atom. 

the  number  of  interstitial  sites  surrounding  a  given  C-site  which  are 
excluded  from  occupancy  (by  other  C  atoms)  if  it  is  filled  by  a  C  atom. 

the  excess  free  energy  of  mixing  of  the  ft  interstitial  solid  solution 
which  is  a  function  of  composition. 
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